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1.  INTRODUCTION 


1 . 1 Object  and  Scope 

The  object  of  this  study  is  to  determine  accurately  the  stresses  and 
displacements  in  two-dimensional  elastic  bodies  with  high  stress  concen- 
trations using  a form  of  the  boundary  integral  equation  technique.  Major 
emphasis  will  be  placed  on  treating  boundary  value  problems  for  bodies 
where  at  some  points  the  radius  of  curvature  of  the  boundary  becomes  very 
small.  Example  of  such  configurations  are  notches  where  the  root  of  the 
notch  is  curved  and  parallel-sided  cracks  of  finite  but  small  width  (wide 
crack).  An  improvement  of  a previously  developed  procedure  for  treating 
the  case  of  an  infinitesimal  crack  provides  a considerable  advantage  in  the 
solution  of  the  wide  crack  problem  as  well. 

The  material  of  the  body  is  assumed  to  be  homogeneous,  isotropic  and 
linearly  elastic.  The  displacements  are  taken  small  enough  that  the  ''inear 
theory  of  elasticity  is  applicable. 

The  technical  importance  of  the  problems  under  consideration  needs  no 
emphasis.  The  purpose  of  the  sample  problems  solved  is  to  display  the 
versatility  and  accuracy  of  the  procedure  rather  than  to  give  numerical 
solution  to  actual  problems  of  engineering  significance. 

1.2  General  Remarks  on  Boundary  Integral  Equation  Method 

The  boundary  integral  equation  (BIE)  method  provides  a powerful  tool 
for  solving  boundary  value  problems  in  potential  theory  and  classical 
elasticity.  Historically,  according  to  Love  [l]1,  the  integral  equation 
method  was  first  applied  in  theory  of  elasticity  by  Betti  in  1872. 

^Numbers  in  brackets  refer  to  entries  in  the  list  of  references. 
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Further  work  was  done  in  this  area  by  Somigliana  in  the  1880's.  The  general 

elasticity  problem1  is  solved  in  the  Somigliana  approach  [2]  by  applying 

Betti's  law  to  the  actual  state  of  stress  and  the  state  of  stress  arising 

2 

from  the  Kelvin's  fundamental  solutions.  The  method  overcomes  the  inherent 
difficulty  of  Betti's  approach  which  requires  construction  of  Green's 
type  functions  for  the  solution  of  the  same  problem. 

The  advent  of  digital  computers  has  been  responsible  for  the  develop- 
ment of  interest  in  the  numerical  treatment  of  elasticity  problems  using 
the  boundary  integral  equation  approach.  The  method  has  been  modified  and 
extended  in  a number  of  investigations  to  solve  two-  and  three-dimensional 
problems  in  isotropic  elastic  media,  and  even  in  elasto-plastic  and  anisotropic 
media  [5-12]. 

The  method  has  also  been  extended  to  the  solution  of  bodies  with 
infinite  stress  gradients  such  as  cracks  and  sharp  re-entrant  corners 
(for  example,  see  Ref.  12-17).  One  method  of  special  interest  was 
developed  by  Barone  and  Robinson  [18].  This  method  makes  use  of  Williams' 
solutions  [19,  20]  with  unknown  coefficients  to  represent  the  stress  and 
displacement  field  in  the  vicinity  of  the  sharp  notches  or  cracks.  This 
work  unlike  other  methods  using  the  Williams'  representation,  employs  new 
auxiliary  solutions  for  the  determination  of  the  unknown  coefficients.  This 
method  will  be  used  extensively  in  this  study.  A more  detailed  explanation 
of  the  method  is  given  in  Section  2.4. 

*In  this  study  we  shall  consider  only  the  so-called  "first"  boundary- value 
problem  in  plane  elasticity  which  is  to  determine  the  displacements  in 
the  interior  of  an  elastic  body  where  the  body  force  and  boundary  tractions 
are  prescribed.  Extension  to  other  well-set  boundary  value  problems  does 
not  appear  to  involve  essential  changes. 

2 

These  solutions  are  sometimes  called  auxiliary  solutions  in  the  BIE  method. 


As  mentioned  before,  the  major  emphasis  will  be  placed  on  treating 


bodies  with  notches  where  the  root  of  the  notch  is  curved  (fillet  problem), 
and  also  bodies  with  a "wide"  crack.  The  solution  of  these  problems 
using  the  ordinary  BIE  method  requires  the  introduction  of  many  points  in 
a small,  localized  region  of  the  boundary.  This  will  cause  the  number  of 
the  unknowns  to  increase  and  contributes  to  poor  conditioning  of  the  final 
equations J Indeed,  in  the  limit,  as  the  configuration  approaches  that 
of  a sharp  notch,  the  method  breaks  down. 

The  overall  stragegy  of  the  proposed  procedure  is  to  modify  the 
process  which  has  been  found  to  work  for  the  sharp  notch.  The  modification 
will  run  into  difficulty  only  if  the  radius  of  the  notch  is  not  small  with 
respect  to  that  of  the  body,  in  which  case  the  ordinary  BIE  method  is 
economically  applicable,  and  one  would  not  speak  of  a stress-concentration 
problem. 

The  method  is  therefore  based  on  the  Barone  and  Robinson  [18]  solution 
technique.  In  that  method  the  Williams'  forms  [19,  20]  are  used  right  up 
to  the  singular  point.  In  the  proposed  modification  a small  body  containing 
the  fillet  is  removed  and  an  extended  Williams'  form  holds  in  an  annular 
region  outside  this  inner  piece  (see  Fig.  4).  The  inner  piece  is  treated 
as  a smooth  body  where  the  details  of  the  configuration  of  the  fillet  are 
taken  into  account  by  the  specified  geometry  of  the  boundary  points.  The 
process  involves  matching  the  smooth  solution  for  the  inner  piece  to  the 

U 

I *It  should  be  noted  that  since  the  equations  of  the  BIE  method  are  not 

banded,  the  computational  penalty,  even  without  regard  to  round-off 
error,  is  severe. 

□ 
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analytic  Williams'  field  on  the  annulus.  The  generalized  displacements 
specifying  the  Williams'  field  are  determined  in  much  the  same  way  as  in 
Ref.  [18],  except  that  the  body  in  question  has  a hole  (the  part  removed 
as  the  inner  body).  After  these  generalized  displacements  are  found,  the 
full  solution  in  the  inner  piece  can  be  constructed  and  the  detailed 
stress  pattern  found  near  the  fillet. 

The  solution  of  the  wide  crack  problem  is  essentially  the  same  as 
for  the  notch  with  a fillet  of  small  radius,  except  that  before  developing 
the  series  solution  for  the  annular  region,  it  is  necessary  to  perturb 
the  Williams'  solution  for  the  corresponding  infinitesimal  crack  in  order 
to  satisfy  the  homogeneous  boundary  conditions  on  the  actual  boundary  of 
the  crack.  Finally,  an  altered  form  of  the  Barone  and  Robinson  solution 
technique  for  the  infinitesimal  crack  is  used  to  avoid  certain  mathematical 
problems,  as  is  explained  in  Appendix  A. 

1 .4  Notation 

The  symbols  used  in  this  study  are  defined  where  they  first  appear. 

For  convenient  reference,  the  more  important  symbols  are  listed  below. 

Some  symbols  are  assigned  more  than  one  meaning,  however  ambiguities  do 
not  occur  in  the  context  of  their  use. 

A , B m-th  eigenfunction 

mm 

Clm*  C2m'  C3m*  m-°’  e19envMtor 

D region  occupied  by  two-dimensional 

body 


.... *4. 
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D 


E 

f 

g 

G(s) 
h,h . 

J 

IL.I(e) 

K 

k 

L 

L 

e 

e,m 

NL 

n(Q) 

P 

P.Q 

i U Q< 

r.r(P.Q) 

i . i(Q> 


D 

I; 

u 


circular  sub-region  of  D with 
radius  e 

displacements  of  perturbed  Williams' 
solution  in  polar  coordinates 

young  modulus  of  elasticity 

density  function 

kernel  of  integral  equations 

= f-g,  product  of  density  function 
and  kernel,  Eq.  (5.1) 

spacing  of  node  points  on  the  boundary 

line  integral 

m-th  generalized  displacement 
boundary  curvature 
boundary  curve  of  region  D 
boundary  curve  of  sub-region  D 
direction  cosine  of  normal  n 
number  of  node  points  in  L 
outward  normal  at  point  Q on  L 
an  interior  point  of  D+L 
boundary  points  on  L 
boundary  point  on  Lc 
distance  between  two  points  P and  Q 
tangent  to  the  boundary  L at  point  Q 

stresses  of  perturbed  Williams' 
solutions 


tj(Q) 

^(P.Q) 

UjW) 

u^lP.Q) 


u , U 

P 4> 


u , U 

p <p 


V \ 

x.y 

z 

z 

€ 

°o 

\n 

p » 4> 


<3  % 0.**  °n< 

( >\)  P ' 


°pp  * ' P$ 


actual  far-fiel d tractions  at 
point  Q 

far-field  traction  at  point  Q due 
to  auxiliary  field  at  point  P 

actual  far-field  displacement  at 
point  Q 

far-field  displacement  at  point  Q 
due  to  auxiliary  field  at  point  P 

displacements  of  the  Williams 
solution  forms 

displacements  of  the  nwdified 
Williams'  solution  or  modified 
perturbed  solution 

net  reciprocal  work  along  L 

rectangular  coordinates 

= x + iy  complex  number 

conjugate  of  z 

radius  of  circular  region  D£ 

included  angle  of  a notch 

m-th  eigenvalue 

polar  coordinates 

stresses  from  Williams'  solutions 
form 

stresses  from  modified  Williams' 
solution  or  modified  perturbed  solution 


The  suiwtion  convention  is  osed  on  the  repeated  subscripts  j in  the 
integral  relationships. 
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2.  THE  BOUNDARY  INTEGRAL  EQUATION  METHOD  AND  ITS  EXTENSION  TO  CORNERS 

2 . 1 Introduction 

As  mentioned  in  Chapter  1,  the  method  of  solution  in  this  study  is 
based  on  an  integral  equation  formulation  of  the  problems.  A brief  review 
of  the  boundary  integral  equation  (BIE)  method  for  bodies  with  smooth 
boundaries  will  be  given  in  Section  2.2.  In  order  to  extend  the  BIE 
method  to  bodies  with  irregular  boundary  point^  it  is  convenient  to  use 
a field  that  represents  the  character  of  displacements  at  that  point. 

Such  an  extension  was  developed  by  Barone  and  Robinson  [18]  and  will  be 
discussed  in  Section  2.4.  Before  the  extension  of  the  BIE  method  to 
irregular  points  can  be  developed,  it  is  necessary  to  know  the  behavior 
of  the  displacements  in  the  vicinity  of  such  points.  The  basic  elasticity 
solutions  which  define  this  behavior,  the  so-called  Williams'  solutions 
[19,20],  are  described  in  Section  2.3.  It  will  be  shown  in  Chapter  3 and 
Chapter  4 how  these  solutions  can  be  modified  to  be  applicable  for  certain 
regions  in  problems  of  notches  with  fillet  of  small  radius  and  in  problems 
with  cracks  of  finite  but  small  width. 

2 • 2 The  Boundary  I ntegral  Equation  Method 

2.2.1  General  Remarks 

In  recent  years  much  attention  has  been  given  to  the  numerical  solution 
of  boundary  value  problems  using  integral  equations.  The  method  has  been 


Irregular  boundary  points  are  those  at  which  the  tangent  to  the  boundary 
does  not  turn  smoothly  or  a sudden  change  in  the  nature  of  the  boundary 
conditions  occurs.  These  include  comers,  notch  tips,  crack  tips,  and  an 
edge  of  a contact  region.  This  is  not  to  be  confused  with  Kellogg's  defini- 
tion of  regular  boundaries  and  surfaces  [37]. 
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used  successfully  for  solving  a wide  range  of  problems  in  applied  mechanics 
(for  example,  see  Refs.  5-23).  The  main  advantage  of  the  boundary  integral 
equation  is  that  for  the  solution  of  a given  problem,  the  unknown  boundary 
values  are  expressed  entirely  in  terms  of  boundary  data.  This  is  in 
contrast  to  the  finite  difference  and  finite  element  methods  which  use 
Interior  field  values  in  addition  to  the  boundary  values  as  unknowns  for  the 
solution  of  the  same  problem. 

The  formulation  of  the  boundary  integral  equation  method  for  the  theory 
of  elasticity  is  based  on  Betti's  law  [l  ].  The  two  sets  of  displacement 
and  corresponding  stress  fields  used  in  Betti's  law  will  be  the  actual 
system  and  another,  auxiliary  system.  The  object  is  determination  of 
unknown  displacements  and  stresses  of  the  actual  system  using  the  known 
fields  of  the  auxiliary  system  as  well  as  whatever  boundary  values  of  the 
actual  system  are  known. 

The  success  of  the  BIE  method  depends  on  the  choice  of  the  auxiliary 
solution.  It  is  desirable  to  have  auxiliary  fields  corresponding  to 
loads  which  pick  out  a component  of  displacement  at  an  interior  point. 

Such  auxiliary  fields  are  those  resulting  from  the  point  loads  acting 
on  a body  in  the  direction  of  the  desired  displacements.  A proper 
choice  of  other  auxiliary  loadings  is  the  key  to  the  convenient  extension 
of  the  BIE  method  to  be  discussed  in  Section  2.4. 

2.2.2  Kelvin  Solution  for  the  Infinite  PI  ate 

Concentrated  force  solutions,  the  so-called  Kelvin  solutions,  are 
auxiliaries  that  will  be  used  extensively  in  this  study.  The  tractions 
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* * 1 

and  displacement  t.j(P,Q)  and  u^(P,Q)  at  a boundary  point  Q due  to  a 
concentrated  force  in  the  i-th  coordinate  direction  at  a point  P are: 

a)  Tractions 
t^(P.Q)-A(rx2  ♦ B)  ~ 

t22(P.Q)=A(ry2  ♦ B)^ 

b)  Displacements 
u*, (P.Q)-Ln(r)  + crx2 

u*2(P,Q)=Ln(r)  ♦ cry2 

where  r = r(P,Q)  and  the  subscripts  x,  y,  n,  s indicate  differentiation 
with  respect  to  the  x,  y,  n,  and  s coordinates  (see  Fig.  1).  The 
constants  A,  B,  and  C for  the  case  of  plane  strain  are: 

Ax(T+vf(3-TvT  ; B4(1'2v)  ; Cs'3-Tv  (2‘3) 

In  Eqs.  (2.3),  E is  Young  modulus  of  elasticity  and  v is  Poisson's  ratio. 

2.2.3  Integral  Equation  for  Regular  Boundary  Points 

Consider  a two-dimensional  isotropic  elastic  body  in  plane  strain 
with  boundary  contour  L (see  Fig.  1).  The  Kelvin  field  resulting  from  a 
concentrated  force  applied  at  a point  P of  the  body  in  the  i-th  coordinate 
direction  will  be  used  as  auxiliary.  Since  the  auxiliary  solution  is 
unbounded  at  P,  a small  region  D.  with  boundary  L(  is  removed  from  the  body. 


tl2(P.Q)-A(rxryrn  * Br$)r 
t2,(P.Q)-A(rxr  r - Brs)J 


VlP.QHVy 

^(P.QXV, 


(2.1) 


(2.2) 


*The  first  script,  (i  * 1,2)  refers  to  direction  of  the  load  and  the  second 
(j  - 1,2)  refers  to  x,  y coordinates. 


I ] 

l . 


Betti's  law  can  now  be  written  for  the  remaining  region  in  the  form: 

tj(Qf)  (P,Q.)ds(Qt ) + j t j (Q ) u*j(P ,Q)ds(Q)= 

”l-L 

t c 

uj(Qe)  t^P.OjdsWJ  + j u.(Q)  t*j(P ,Q)ds(Q) 

L-L 


1 


(2.4) 


where  Q represent  the  point  on  the  circular  boundary  L The  kernels 

t*.(P,Q ) and  u..(P,Q)  are  defined  in  the  previous  section.  The  u,(Q) 

1 j i j J 

and  t . ( Q ) are  the  boundary  displacement  and  traction  at  point  Q.  After 
J 

taking  the  limit  c >0.  we  obtain 


°ijuj(p)’ 


tj(Q)  u..(P.Q)ds 


L 

L J 


(Q)  ti.(P,Q)ds. 


(2.5) 


1 im 
e‘0 


lim 

c»0 


/ t.(Qf)  u*j (P *Qe )ds (Qe )=0, 

c 

l t;j(P.Qt)ds(Qc)-«ijuJ(P) 

e 


(2.6) 

(2.7) 


In  Eq.  (2.5)  the  integrals  are  to  be  understood  in  the  Cauchy  principal 
value  sense.  The  displacement  vector  at  an  interior  point  P is  denoted 
by  uj(p).  The  value  of  the  for  an  Interior  point  or  a regular  boundary 
point  is: 


ay  3 0 for  i f j 

c*ij  » - 0g  A(B  + 0.25)  for  i « j 


(2.8) 


^n  the  equations  of  this  section  the  sunmation  convention  is  used  for 
lower  case  subscripts. 
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where  oQ  * 2*  for  an  Interior  point  and  «g  " " for  a point  on  a smooth 
boundary.  The  value  of  for  a sharp  corner  Is:* 

a1j  ■*  (A/4)[cos(2^q  2eQ)-  cos(2bQ)]  for  1 t j (2  g) 

'U  --{tf0A(B  + 0.25)  ♦ (A/4)[  s1n(2t'0  + 2eQ)-  s1n(2BQ)]1  for  1 - j 

where  the  positive  sign  holds  for  j*l  and  negative  sign  for  j«2.  The  angle 
i*q  is  the  Internal  angle  of  the  notch  and  t'g  Is  shown  In  Fig.  2. 

2 . 3 The  Willi ams ' So luti ons 
2.3.1  I n t rod uction 

In  the  development  of  the  proposed  method  of  this  study,  extensive 
use  will  be  made  of  the  displacement  field  in  the  neighborhood  of  an 
irregular  boundary  point.  Analytical  results  for  determination  of  the 
character  of  the  displacement  and  stress  fields  are  well  known  [19,20]. 

A more  recent  approach  by  0.  K.  Aksentian  [24],  which  is  based  on  a 
three-dimensional  elastic  analysis  of  the  problem,  provides  a somewhat 
more  general  solution  of  the  problem  in  the  form  of  asymptotic  expressions. 

A brief  description  of  the  Williams'  solution  for  two-dimensional  bodies 
In  plane  strain  will  be  given  hero. 

The  displacement  field  in  the  neighborhood  of  an  irregular  point  in 
polar  coordinates  may  be  written  in  the  form 

“o  ' AU)'’'  (2.10) 

u,  • BU)/ 

Hhe  Kelvin  auxiliaries  are  used  to  pick  out  only  the  rigid  body  displacements 
at  a sharp  corner  (see  Section  3.5). 


-* 
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Solving  the  differential  equations  of  equilibrium  for  the  functions 
) and  B($)  gi ves 

A(*)  = (x-3+4v)[c^cos(x-l  )4>+c2sin(x-l  )4>  ]+c3cos  (x+1  H+c^sinU+l  )$ 

(2.11) 

B(^)=(x+3-4v)[-c1sin(x-l  H+c2cos(x-l  )$]-c3sin(x+l  H+c4cos(x+l  H 

where  the  c's  are  arbitrary  constants  of  integration  and  vis  Poisson's  ratio. 
The  corresponding  stresses  in  polar  coordinates  are: 

o ■ PU)pX_1 

pp 

= (2.12) 

o 3 R(  4>  )p X 1 
p<fr 

where 


P (♦  )=2uX{  (x-1  )[c1cos(x-l  )^+c2sin(x-l  )$]+c3cos(x+l  H*c4sin(x+1  )*} 

QU  )”~2uX{  (x+l)[Cj cos  (x-1  )<j>+c2sin(x-l  )«|>]+c3cos (x+1  )<j>+c4sin(x+l  )<j>)  (2.13) 

R(4>)32pX{  (x-l)[-c^sin(x-l  )$+c2cos(x-l)$]-c3sin(x+l  )<i>+c4cos(x+l  )<j>> 

and  y is  the  shear  modulus. 


The  constants  c^  must  be  determined  so  that  displacements  and  stresses 
satisfy  the  imposed  boundary  conditions.  If  the  body  is  a wedge,  boundary 
conditions  will  be  given  along  the  sides  of  the  wedge.  For  example, 
specifications  of  two  homogeneous  tractions  boundary  conditions  along 
each  edge  leads  to  an  eigenvalue  problem  [19,24]  with  characteristic 

equation  of  the  form 

2 2 2 
sin  xe^x  sin 


(2.14) 
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where  eQ  is  the  internal  angle  of  the  wedge.  There  are  an  infinite  number 
of  eigenvalues  ( = 1.2, 3,..)  which  satisfy  Eq.  (2.14)  and  lead  to  an 
infinite  number  of  eigenfunctions,  each  of  which  is  determined  to  within 
a multiplicative  constant. 

2.3.2  Series  Representation  of  the  Displacement  Field 

It  was  shown  that  there  are  an  infinite  number  of  displacement  fields 
each  of  which  satisfies  the  homogeneous  boundary  conditions.  In  general,  the 
results  occur  in  complex  conjugate  pairs  of  solutions,  each  pair  leading  to 
two  real  solutions.  Detailed  representation  of  these  fields  are  given  by 
Barone  [18].  For  simplicity  of  the  arguments,  the  eigenfunctions  will  be 
considered  in  the  form  p%(<*>)  even  if  \ is  complex. 

Any  combination  of  the  Williams'  solutions  is  also  a solution.  That 
is,  the  following  expressions  are  possible  solutions  of  the  homogeneous 
problem. 


u 

p 


= l 


K u 
m pm 


= l 


K u 
m $m 


o 

pp 


l K 


m 'ppm 


- I K 

m 


m°<f><f>m 


o 


P 


1 K 


map<j>m 


u m = A ( (j>)PArn 
pm  m * M 

uAm  = Bm(4>)pAm 
4>m  m 

a = P ( 4>  )p Afn  ^ 
ppm  m 


(2.15) 
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In  the  Barone  and  Robinson  solution  techniques  [18],  the  K's  serve 
as  generalized  displacements.  In  order  to  insure  that  the  displacements  in 
Eqs.  (2.15)  are  bounded,  the  real  part  of  each  eigenvalue  xm  must  be 
greater  or  equal  to  zero.  Such  solutions  will  be  termed  "proper" 

Williams'  solutions.  However,  in  Eq.  (2.14)  if  xm  is  a solution,  then 
-\  is  also  a solution.^  A solution  with  the  real  part  of  the  eigen- 

2 

value  less  than  zero  will  be  called  an  "improper"  Williams'  solution. 

As  will  be  seen,  these  improper  solutions  are  used  as  auxiliaries  for 

the  determination  of  the  generalized  displacements.  They  can  be  given  as 

* 


= I *V 

£ m pm 
m 

* 

★ 

U m = 
pm 

AmU)pXm 

★ 

= l 

£ m <j>m 
m 

* 

v- 

B*U)pXnf] 

★ 

_ * * 

= y k o 

jj,  ro  Wro 

. 

★ 

o 

= P ( 4>  )p  Xm 

ppm 

m 

★ 

„ * * 

= T K o 
' ^ m <j>4>m 

* 

* 

a 4>  4>m 

= Q*(<t>)pX,n"1 

= i kV 

_ m p 

m 

9 

* 

ap4>m 

= ^(♦)pA"“1 

(2.16) 


The  proper  Williams'  solutions  appear  to  provide  a complete  set  for 
expanding  an  arbitrary  distribution  of  self-equilibrated  tractions  on  an 
artificial  circular  boundary  in  a body,  having  the  irregular  point  as  its 
center.  The  numerical  verification  of  this  has  been  carried  out  in  several 


It  is  shown  in  Ref.  18  that  solution  exists  for  these  eigenvalues  which 
satisfy  the  homogeneous  boundary  conditions  for  well-set  problems. 

'Improper  Williams'  solutions  correspond  to  unbounded  strain  energy  near 
the  singular  point. 
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cases.  For  example,  let  us  take  a disk  with  a crack  such  that  tip  of  the 
crack  is  at  center  of  the  disk  (see  Fig.  16).  For  an  arbitrary  self- 
equilibrated  tractions  on  the  boundary  of  the  disk  one  can  take  a large 
enough  number  of  eigenfunctions  in  an  expansion  to  reduce  the  error  on 
the  boundary  of  the  disk  below  any  prescribed  value.  As  a numerical 
example,  the  disk  is  loaded  by  radial  tension  and  Eqs.  (2.15)  were  used 
for  the  expansion.  Table  1 gives  the  coefficients  of  the  expansion  and 
Table  2 the  computed  tractions  for  points  on  the  boundary.  The  results 
will  be  discussed  fully  in  Chapter  5. 

It  seems  that  the  circular  shape  of  the  boundary  is  not  an  unimportant 
factor  in  the  convergence  of  the  solution,  for  the  results  are  not  satis- 
factory on  some  bodies  with  non-circular  boundaries,  using  say,  a collocation 
method  [25].  It  is  also  shown  by  Barone  and  Robinson  [18]  that  the  displace- 
ment field  of  the  irregular  point  is  valid  in  a small  region  about  the 
irregular  point.  One  might  expand  the  region  by  adding  to  the  number  of 
eigenfunctions  in  the  expansion  expression  and  still  retain  accuracy. 

However,  the  region  for  which  the  expansion  is  valid  can  not  include  other 
singular  points  on  the  boundary  of  the  body.  These  are  inevitable  at  any 
corner  or  any  point  where  the  boundary  condition  changes  suddenly.  These 
observations  make  it  reasonable  to  limit  the  area  where  the  Williams'  series 
is  used  to  a fairly  small  region  around  its  generating  singularity. 

In  general,  the  proper  and  improper  Williams'  solution  taken  together 
(the  "extended"  Williams'  solution)  will  be  needed  to  represent  a field  in 
an  annulus  having  two  circular  boundaries.  The  analogy  with  spherical  and 
cylindrical  harmonics  suggests  that  singular  solutions  must  be  considered 
when  the  origin  is  not  part  of  the  body. 

D 
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2 • 4 The  Solution  Technique  of  Barone  and  Robinson 


The  formulation  of  the  boundary  integral  equations  for  elastic  bodies 
with  irregular  boundary  points  is  given  by  Barone  and  Robinson  [18].  In 
this  method  the  proper  Williams'  solution  Eqs.  (2.15),  with  K's  as  unknown 
generalized  displacements  are  used  to  represent  the  displacement  field  ir 
the  vicinity  of  the  irregular  boundary  point  up  to  so-called  transition 
points.  The  integral  equation  for  the  regular  boundary  points  are  then 
written  as  described  in  Section  2.2.3.  For  irregular  boundary  points 
a new  kernel  was  derived  which  turns  out  to  be  the  eigenfunctions  of  the 
improper  Wi 11 iams 1 solution  Eqs.  (2.16).  These  auxiliary  solutions  are  used 
to  pick  out  the  unknown  generalized  displacements.  However,  for  the  first 
two  K's  corresponding  to  rigid-body  displacements  the  Kelvin  auxiliary 
solutions  are  used. 

Since  the  improper  Williams'  solutions  are  singular  at  an  irregular 
point,  a small  circular  piece  of  radius  e containing  the  irregular  point 
removed  from  the  body  before  the  reciprocal  theorem  is  applied  to  the  actual 
and  auxiliary  solution.  The  reciprocity  relation  for  the  body  after  taking 
the  limit  e -*  0 is 


M 

l 

m=3 


- K + I t..(P,Q)  u . (Q)  ds(Q)  = 
im  m /,  ij  j 


(2.17) 


u*.(P,Q)  tj (Q)  ds(Q) 


where  in  Eq.  (2.17)  the  standard  summation  convention  is  used  on  the  lower 
case  subscripts.  The  t*.(P,Q)  and  u..(P,Q)  are  tractions  and  displacements 


. ■ 


of  the  i-th  improper  Williams'  solution.  The  t . (Q ) and  u - (Q ) are  the  traction 

J J 

and  the  displacement  of  the  boundary  point  Q in  the  j-th  coordinate  direction 
and  P is  the  irregular  point  at  which  the  auxiliary  solutions  are  singular. 
Here  M is  the  total  number  of  proper  Williams'  solutions  used,  including 
the  three  rigid-body  motions.^  Furthermore, 

a.  = 0 for  i i m 

im 

a.  = F.  for  i = m , m > 3 

im  im 

where  F^m  is  given  by  Eq.  (C.5). 

As  can  be  seen,  the  m-th  auxiliary  solution  pick  out  the  m-th 
generalized  displacement  and  relates  it  to  the  far  field  boundary  values. 

The  ordinary  integral  equations  at  regular  points  along  with  these  new  ones 
form  a coupled  system  involving  the  generalized  displacements  and  other 
boundary  values  as  unknowns.  This  leads  to  a strongly  diagonal  system  of 
linear  algebraic  equations  which  is  solved  for  the  unknown  values.  As 
mentioned  before,  the  generalized  displacements  depend  on  the  boundary 
values  at  fairly  remote  points.  This  appears  to  contribute  to  a numerically 
stable  solution.  Physically,  the  generalized  displacements  must  depend 
on  the  applied  loads  in  the  far  field.  The  fact  that  this  dependence 
is  explicit  in  the  numerical  procedure  is  one  advantage  of  the  method. 


Hhe  two  translations  are  "picked  out"  by  Kelvin  solutions.  The  rigid- 
body  rotation  is  in  the  set  of  Williams'  solutions. 
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3.  APPROXIMATE  DETERMINATION  OF  STRESSES  AND 
DISPLACEMENTS  NEAR  A NOTCH  WITH  A SMALL  FILLET 

3.1  I nt reduction 

The  main  objective  of  this  Chapter  is  to  develop  an  effective  numerical 
scheme  for  the  determination  of  stresses  and  displacements  at  points  in  a 
small  region  containing  the  fillet  (the  inner  region).  This  is  done  by 
developing  a set  of  solutions  for  the  inner  region  and  tying  them  to  the 
rest  of  the  body  by  expanding  the  solutions  in  a series  with  unknown 
coefficients.  This  procedure  allows  us  to  overcome  the  poor  conditioning 
of  equations  which  would  arise  by  including  the  details  of  the  boundary  of 
the  small  fillet  in  the  same  calculation  as  the  boundary  of  the  large  body. 

As  will  be  seen,  the  proper  Williams ' solutions  (Eqs.  2.15)  corres- 
ponding to  the  sharp  notch  are  sufficient  for  developing  the  inner  solutions, 
but  the  extended  (full)  Williams'  solutions  must  be  used  for  the  series 
expansion  in  the  so-called  annular  region. 

3 . 2 The  Solution  of  the  Inner  and  Annular  Regions 

We  take  a two-dimensional  body  with  a notch  where  the  root  of  the  notch 
is  curved  (see  Fig.  3).  We  shall  consider  three  regions  of  the  body: 

i.  An  inner  piece  containing  the  fillet  and  bounded  by  a circle 
centered  at  the  Intersection  of  the  two  straight  edges  of  the 
notch  (Point  P,  Fig.  4); 

ii.  An  annular  region  adjacent  to  the  inner  piece  having  as  its 

outer  boundary  a circle  passing  through  the  transition  points;  and 
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Hi.  The  outer  body,  or  the  large  body,  which  consists  of  the  rest  of 
the  body  exterior  to  the  annulus  and  the  annulus  itself. 

Let  us  analyze  the  inner  body  by  the  ordinary  BIE  method  as  described 
in  Section  2. 2. 3, for  one  of  the  proper  Williams'  tractions  applied  to  Its 
outer  boundary.  A field  which  is  determined  In  this  way  for  the  inner 
body  and  its  boundary  Is  termed  a "unit  field."  Now  the  displacements 
as  well  as  the  tractions  are  available  on  the  circle  ABC  (see  Fig.  4) 
separating  the  inner  body  from  the  annulus.  As  analytic  solution  consist- 
ing of  Williams'  solutions  is  sought  in  the  annulus  which  matches  the 
displacements  and  stresses  on  ABC  and  thus  completely  guarantees  continuity 
of  displacements  and  tractions  between  the  inner  body  and  the  annulus.  In 
general  this  will  require  improper  as  well  as  proper  Williams'  solutions. 

To  this  end,  let  us  expand  the  computed  field  for  the  annular  region 
(ABC-FED,  Fig.  5)  in  a series  using  the  proper  and  Improper  Williams' 
solutions  (Eqs.  2.15  and  2.16),  which  were  heuristical ly  justified  in 
Section  2.3.  What  results  is  termed  a "modified  Williams'  solution."  It 
is  modified  by  the  presence  of  the  fillet,  for  if  the  notch  were  sharp, 
the  result  would  have  been  the  single  proper  Williams'  solution  corresponding 
to  the  applied  tractions. 

It  might  seem  incorrect  to  use  only  proper  Williams'  tractions  to 
find  the  unit  fields  and  then  to  match  these  unit  fields  to  the  field  on 
the  annulus  with  both  improper  and  proper  solutions.  However,  as  was 
mentioned  in  Section  2.3,  any  traction  on  the  circle  can  be  expanded  in 
proper  Williams'  tractions.  The  Inclusion  of  the  improper  ones  in  addition 
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would  lead  to  non-uniqueness  of  the  coefficients.  It  should  be  remembered 
that  the  series  representation  in  the  annulus  matches  displacements  as 
well  as  tractions. 

The  general  solution  for  the  annular  region  is  now  expressed  in  a 
series  of  the  modified  Williams'  solutions,  the  unknown  coefficients  of 
which  are  generalized  coordinates  for  the  outer  body.1  Once  these  coeffi- 
cients are  determined,  the  solution  for  the  inner  body  is  just  a linear 
combination  of  the  previously  calculated  inner  solutions  (unit  fields). 

3. 3 Calculation  of  the  Modified  Williams'  Solution 

3.3.1  Series  Representation  of  the  Solution 

As  described  earlier,  in  order  to  develop  a new  solution  to  represent 
the  displacement  field  of  the  annular  region  it  is  necessary  to  remove  the 
inner  piece  from  the  body.  This  piece  is  loaded  by  tractions  corresponding 
to  one  of  the  proper  eigenfunctions  of  the  Williams'  solution  for  the 
sharp  notch,  say  the  m-th  one.  The  pointwise  displacements  are  then 
computed  using  the  ordinary  BIE  method.  This  process  is  carried  out  for 
the  first  M eigenfunctions  of  the  Williams'  solution.  Therefore,  corre- 
sponding to  each  eigenfunction  a new  field  is  developed  which  reflects  the 

*At  first  sight,  the  procedure  would  seem  to  be  open  to  Hadamard's  objection 
to  solution  of  an  elliptic  problem  with  Cauchy  initial  data  [29,  32].  It 
should  be  remembered,  however,  that  here  the  coefficients  are  finally 
determined  by  the  conditions  at  the  outer  boundary  of  the  body.  All  that 
has  been  done  is  to  partition  the  total  extended  set  of  the  Williams' 
solutions  in  a way  convenient  for  calculation.  It  is,  however,  important 
that  the  number  of  proper  Williams'  solutions  (the  ones  that  grow  with  p) 
be  equal  to  the  number  of  K's.  This  is  automatic  in  the  present  procedure. 
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Influence  of  the  fillet  on  the  displacement  and  stress  fields.  The  results 
of  each  solution  are  now  represented  In  series  of  the  form: 
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where  the  a^'s  are  the  unknown  coefficients  and  the  expansion  functions 
are  the  eigenfunctions  corresponding  to  the  sharp  notch.  They  include 
the  real  as  well  as  the  complex  eigenfunctions.  However,  here  the  bar  on 
the  displacements  or  stresses  does  not  represent  the  complex  conjugate, 
rather  the  notation  is  used  to  represent  the  modified  Williams'  solution 
for  the  annular  region.  The  scripts  j<0  are  used  to  represent  the  eigen- 
functions with  the  eigenvalues  X 0;  j>0  means  Xj>0.  The  relation  |j|>4 
indicates  that  the  three  rigid  body  motions  are  excluded  in  the  expansion. 
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3.3.2  Matching  of  the  Field  of  the  Inner  and  Annular  Regions 

Betti's  law  is  used  for  the  determination  of  the  a^'s,  the  unknown 
coefficients  of  the  m-th  modified  Williams'  solution.  The  objective  is 
to  connect  the  fields  of  the  inner  and  annular  regions  smoothly. 

Consider  a smaller  annulus  in  the  body  with  outer  radius  R2 
(ABCA'B'C',  Fig.  5).  This  piece  is  loaded  along  the  boundary  ABC  with 
the  tractions  of  the  m-th  unit  field  of  the  inner  piece,  and  along  the 
boundary  A'B'C'  with  the  tractions  of  the  series  representation  of  the  same 
solution.  The  k-th  improper  Williams'  solution  is  now  chosen  as  an 
auxil iary. 

Betti's  law  applied  to  this  annular  piece  results  in 
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where  Rq  is  the  radius  the  inner  piece  and  the  Fmk  is  given  by  Eq.  (C.5) 
in  Appendix  C.  It  is  also  shown  in  that  Appendix  that: 


Fjk  = 0 for  j f k or  xk  + xj  ^ 0 

* * 

Fjk  t 0 for  j = k or  xk  + Xj  = ® 
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Therefore, 
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As  can  be  seen,  the  a^'s  are  functions  of  the  outer  radius  RQ,  of 
the  inner  piece  only.  A question  naturally  arises:  What  should  be  the 
size  of  the  Inner  piece?  As  mentioned  earlier,  the  Inner  piece  must  be 
large  enough  to  Include  the  whole  fillet  so  that  the  boundaries  In  the 
annular  region  are  straight.  This  will  Insure  homogeneous  boundary  condi- 
tions for  the  edges  of  the  annular  region.  The  Inner  piece  should  also  be 
small  enough  that  the  stress  gradients  are  not  much  larger  than  the  maximum 
stress  divided  by  the  diameter.  This  should  guarantee  accurate  results  by 
the  ordinary  numerical  Integral  equation  process. 

In  this  stucty  the  radius  of  the  Inner  piece  Is  taken  to  be  2.5  times 
of  the  radius  of  the  fillet. 

3.4  Series  Representation  of  the  Displacement  Field  for  the  Annulus 

Once  the  unknown  coefficients  of  the  modified  Williams1  solution  are 
determined,  we  have  a solution  for  the  annular  region  which  takes  account 
of  the  fillet  of  the  Inner  piece.  This  solution  can  be  written  In  series 
form 
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where  the  K's  are  generalized  displacements  and  the  expansion  functions 
are  the  modified  Williams'  solution.  Here  M Is  the  total  number  of  modified 
functions  for  the  solution  of  a specific  problem.  Note  that  the  first  three 
functions  In  the  series  are  the  displacement  field  for  rigid  body  motions 
which  are  Identical  for  both  modified  and  original  Williams'  solutions. 

3.5  Solution  for  the  Large  Body  (Loaded  Body) 

The  solution  procedure  for  the  large  body  Is  essentially  the  same  as 
was  described  in  Section  2.4  for  the  solution  of  the  body  with  a sharp 
notch  using  the  BIE  method.  In  fact,  the  same  Improper  Williams'  solutions 
used  as  auxiliaries  for  the  sharp  notch  problem  also  serve  as  auxiliaries 
for  the  large  body  formed  by  the  removal  of  the  Inner  piece J These 
auxiliaries  will  pick  out  the  unknown  generalized  displacements  given 
In  Eqs.  (3.7). 

The  Inner  piece  may  be  considered  as  the  counterpart  of  the  piece  of 
radius  e centered  at  the  Irregular  point  in  the  sharp  notch  problem.  This 
infinitesimally  small  piece  had  to  be  taken  out  from  the  body  to  insure 
applicability  of  Betti's  law. 

For  evaluation  of  integrals  along  the  boundary  formed  by  the  removal  of 
the  finite  Inner  piece  (Fig.  6),  an  equivalent  path  of  Integration  can  be 


^ere,  just  as  the  case  of  the  sharp  notch,  the  Kelvin  auxiliaries  are 
used  to  pick  out  the  rigid  body  translations. 
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chosen  In  order  to  simplify  the  computation  of  the  Integrals.  One  simple 
path  of  integration  consists  of  the  boundary  edges  of  the  fictitious  sharp 
notch  formed  by  the  extension  of  the  two  straight  parts  of  the  boundary  of 
the  inner  piece  (see  Fig.  7).  The  following  consideration  will  show  that 
this  is  an  equivalent  path  of  integration.  Of  course.  It  should  be  observed 
that  for  the  actual  body  the  modified  Williams'  solutions  developed  in 
Section  3.4  are  not  valid  for  any  point  inside  the  inner  body.  However, 
these  solutions  do  satisfy  the  equilibrium  equations  and  are  continuous 
in  the  whole  region  except  possibly  at  the  tip  of  the  fictitious  notch. 

Now  let  us  write  Betti's  work  integral  for  the  body  ABCDEF  of  Fig.  7, 
using  as  the  applied  field  the  modified  Williams'  solutions  and  as  the 
auxiliary  field  the  displacements  and  tractions  from  a source  outside 
the  body.  We  then  have 


!ABC  + rCDEFA  = 0 ; *ABC  "_ICDE FA 


(3.8) 


The  boundaries  ABC  and  CDEFA  are  shown  in  Fig.  7 and 


[Vj  - ujV* 


(3.9) 


Since  some  of  the  solutions  are  unbounded  at  point  E a small  piece 
with  radius  e centered  at  this  point  is  removed  from  the  body  before 
application  of  the  Betti's  law  to  the  larger  body.  Details  of  the 
approximation  of  the  integral  expressions  and  determination  of  the  unknowns 
are  given  in  Chapter  5. 
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Once  the  Integral  equations  are  approximated  and  the  final  linear 
equations  solved  for  the  unknown  generalized  displacements  K's  and  other 
unknown  boundary  values,  the  displacements  and  stresses  at  any  point  In 
the  annulus  can  be  determined  using  Eqs.  (3.5).  The  unknown  fields  for 
points  In  the  Inner  piece  can  be  determined  by  applying  Eq.  2.5  with 
already  computed  boundary  data. 
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4.  APPROXIMATE  DETERMINATION  OF  STRESSES  AND 
DISPLACEMENTS  NEAR  THE  TIP  OF  A PARALLEL-SIDED 
CRACK  OF  FINITE  WIDTH 

4. 1 Introduction 

The  main  object  of  this  chapter  is  to  develop  a numerical  scheme  for 
evaluation  of  stresses  and  displacements  at  points  close  to  the  tip  of  a 
parallel-sided  crack  of  finite  but  small  width.  For  the  sake  of  brevity, 
this  problem  will  be  called  the  "wide  crack"  problem.  The  main  features 
of  the  procedure  will  closely  follow  those  described  in  Chapter  3 for  the 
notch  with  a small  fillet. 

There  is,  however,  one  significant  difference  between  the  problem  of 
the  notch  with  small  fillet  and  the  wide  crack  problem.  The  Williams' 
solutions  used  in  Chapter  3 satisfy  the  homogeneous  boundary  conditions 
on  the  notch  boundary  in  the  annular  region.  By  contrast,  here,  the 
Williams'  solutions  corresponding  to  the  infinitesimal  crack  give  non-zero 
tractions  on  the  edges  of  the  wide  crack.  In  order  to  satisfy  the  homogen- 
eous boundary  conditions  along"  the  edges  of  the  wide  crack,  the  tractions 
will  be  removed  by  a boundary  perturbation  scheme.  What  results  will  be 
termed  a "perturbed  Williams'  solution"  for  the  wide  crack.  These  solutions 
will  play  the  same  role  in  the  wide  crack  problem  as  the  ordinary  Williams' 
solutions  played  in  the  problem  of  the  notch  with  a small  fillet. 
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1 

4.2  Determination  of  the  Perturbed  Williams1  Solutions 
4.2.1  Boundary  Conditions  for  the  Perturbation  Functions 

An  ordinary  or  unperturbed  Williams'  solution  leaves  tractions  on  the 
actual  boundary  because  the  traction  boundary  conditions  were  imposed  at 
the  center  line  of  the  crack,  not  on  the  actual  crack  edges.  A perturbation 
sequence  of  solutions  will  be  introduced  to  remove  these  undesired  tractions 
on  the  real  crack  edges.  The  perturbation  functions  will  be  found  as  the 
solution  of  boundary  value  problems  with  the  center  line  of  the  crack  as 
the  boundary,  in  much  the  same  way  as  the  homogeneous  boundary  tractions 
for  the  unperturbed  Williams'  solutions  were  imposed  on  the  center  line 
of  the  crack.  It  should  also  be  mentioned  that  the  unperturbed  Williams' 
solutions  are  eigenfunctions  (see  Section  2.3),  with  respect  to  these 
center-line  boundary  conditions.  Some  care  must  then  be  exercised  in 
avoiding  those  perturbation  functions  which  only  change  the  magnitude  of 
the  any  eigenfunction  present  in  the  entire  solution.  This  will  keep  the 
perturbations  as  small  as  possible. 

For  clarity  and  computational  ease,  symmetric  and  anti -symmetric 
cases  are  treated  separately.  This  will  make  it  possible  to  deal  with 
just  a single  boundary.  Perturbation  of  the  Williams'  solutions  will  be 
such  that  the  synmetric  or  anti -symmetric  character  of  the  original  solution 
is  preserved.  I 


It  is  convenient  to  use  the  single  complex  function  * = S + iS 

yy  y 

to  represent  two  components  of  the  stress  field  of  the  perturbed  solution. 
The  relations  $(x,i  d)  = 0 for  x < 0 express  the  desired  boundary  conditions 
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for  the  wide  crack  problem.  Let  us  write  the  function  » In  the  form. 


* - l (4-D 

n*0  n 

where  M Is  total  number  of  perturbations  and  *n  gives  the  stresses 
corresponding  to  the  n-th  perturbation.  Here  6 Is  a dimensionless 
measure  of  width,  equal  to  the  half  width  of  the  crack  d,  divided  by 
the  radius  of  the  Inner  body  (as  in  Section  3.2). 

Let  us  express  the  stresses  along  the  boundary  edge  x < 0 and  y - d 
by  expanding  Eq.  (4.1)  in  a power  series  about  y = 0 for  x < 0. 


*0(x,6)  = 

6 3* 

♦0  +-  — + 
u 1!  3y 

2 2 

6C  3£*c 

2!  17 

a.^x.a)  = 

5.*1  + 

S1  !!i 
1!  3y 

s2.*2(x.«)  * 

62.*2 

3! 


2! 

1! 


(4.2) 


In  Eqs.  (4.2)  all  the  functions  and  their  derivatives  are  evaluated 
along  the  half  line  y = 0,  x < 0.  The  function  corresponds  to  one  of 
the  Williams'  solutions  for  the  infinitesimal  crack.  It  is  either 
symmetric  or  antl-symnetrlc  and  Is  determined  so  that  $ (x,0)  * 0 for 
x < 0.  In  order  to  determine  the  boundary  condltons  for  the  subsequent 
functions  in  the  perturbation,  we  set 
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The  boundary  condition  for  the  function  4^  can  now  be  determined  using 
the  first  equation.  After  satisfying  the  given  boundary  condition  is 
found,  the  second  equation  can  be  used  to  determine  the  boundary  conditions 
for  the  4>2»  and  so  on. 

It  should  be  not’d  that  the  procedure  as  described  defines  the  boundary 
conditions  for  the  perturbation  functions  on  the  edges  of  the  corresponding 
infinitesimal  crack,  rather  than  on  the  actual  boundary.  This  will  permit 
determination  of  the  perturbation  functions  in  forms  similar  to  those 
encountered  in  Chapter  2 in  the  solution  of  the  Williams'  problem  for  an 
infinitesimal  crack.  The  main  difference  is  that  for  the  perturbation 
functions,  solutions  will  be  sought  with  specified  non-zero  values  on  the 
half  line  4>  = ±tt  as  opposed  to  the  zero  tractions  of  the  Williams'  case. 

4.2.2  Determination  of  the  Perturbation  Functions  from  their  Boundary 
Conditions 

The  $n  functions  (Eq.  4.1)  have  to  be  determined  subject  to  the  boundary 
conditions  given  by  Eqs.  (4.4).  The  general  form  of  the  function  will  be 
taken  as 

‘k  -ji,  Bjk*jk  ' (4-5) 

The  e's  are  unknown  coefficients  and  the  ip ' s are  certain  general 
solution  forms  (see  Appendix  E).  The  general  solution  forms  that  will  be 
used  include  solutions  similar  to  those  of  Williams'  type.  However,  it 
turns  out  that  these  solutions  alone  cannot  satisfy  the  boundary  conditions 
given  by  Eqs.  (4.4),  requiring  the  introduction  of  other  solutions. 
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The  difficulty  may  be  seen  by  the  following  considerations.  A 

perturbation  solution  must  in  general  satisfy  two  real  non-homogeneous 

boundary  conditions  (Eqs.  4.4),  on  the  line  = ±n.  The  first  perturbation 

x -2 

function,  $,  has  boundary  tractions  varying  as  p m . For  the  infinitesimal 

crack  problem  xm is  also  an  eigenvalue  if  xffl  is.  This  means  that  there 

are  tractions  varying  as  p m in  an  eigenfunction  solution  and  as  a 

consequence,  there  is  only  one  non-zero  boundary  value  that  can  be  removed 

using  the  single  non-eigenfunction  solution.^  Another  form  of  solution 

is  then  needed.  Considering  xm  as  a parameter,  the  general  Williams'  form 

(having  no  pre-assigned  boundary  conditions)  can  be  differentiated  with 

2 

respect  to  the  parameter  xm  to  ottain  a new  solution  form  . The  variation 

of  tractions  with  p in  this  new  solution  Involves  a logarithmic  term 
x -2 

p m Ln(p).  Note  that  two  symmetric  and  two  anti-symmetric  solutions  have 

arisen,  not  one  of  each  as  was  wanted.  However,  the  p m Ln(p)  traction 

x -2 

must  be  removed  at  the  same  time  as  the  p m , for  these  are  comparable  in 

magnitude.  It  turns  out  that  this  can  be  done  only  if  the  coefficients 

of  the  logarithmic  term  are  in  the  same  proportion  as  the  eigenfunction 

for  xm-l  (in  the  displacement).  Thus  a single  new  solution  results  which 

x -2 

can,  in  fact,  remove  the  p m traction  and  introduces  no  new  traction 

x 2 

going  like  p m_tLn(p)  on  the  boundary  # = t*. 


We  have  divided  the  Williams'  solutions  corresponding  to  the  infinitesimal 
crack  into  two  groups,  one  symmetric  and  the  other  anti -symmetric.  For 
each  case  we  have  two  constants  of  integration  (the  c's  in  Eqs.  2.11)  in 
our  disposal.  Since  one  of  the  two  independent  solutions  that  can  be  derived 
is  an  eigenfunction,  there  remains  only  one  non-eigenfunction. 


2 


See  Appendix  B. 
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In  the  second  perturbation  four  functions  will  be  needed  to  remove 
the  tractions  going  like  p m and  p m Ln(p)  from  the  boundary  edge  = n. 

Two  of  the  functions  will  be  similar  in  form  to  those  used  in  the  first 
perturbation  (with  an  exponent  less  by  one).  The  other  two  functions  have 
to  remove  the  boundary  tractions  going  like  pXm'\n(p).  A second  logarithmic 
solution  will  also  be  used;  here  the  coefficients  of  the  logarithmic 
part  will  not  he  taken  as  an  eigenfunction.  For  the  fourth  solution 
form,  a second  differentiation  of  the  general  Williams’  form  with  respect 
to  is  needed.  This  solution  involve  terms  going  like  pXnr3,pXm~\n(p ) and 
p m Ln  (p).  Just  as  in  the  case  of  the  first  perturbation,  the  coefficients 
of  the  last  term  are  taken  to  be  proportional  to  the  eigefunction  corres- 
ponding to  xm-2  (in  displacement)1;  the  process  continues  in  this  fashion. 

In  general  for  the  n-th  perturbation,  2n  functions  will  be  needed  in  order 
to  remove  the  tractions  on  the  boundary  going  like  Lnr(p),(r  = 0,  1,  2...,  n-1). 

As  an  example,  the  solution  forms  and  their  participation  coefficients 
for  the  first  perturbation  are  given  in  Appendix  E. 

4.3  Determination  of  the  Modified  Perturbed  Solutions 

Once  the  perturbed  Williams’  solutions  are  determined,  a small  inner 
piece  of  the  body  surrounding  the  crack  tip  (see  Fig.  9)  is  removed  from 
the  body.  This  piece  is  loaded  by  the  tractions  corresponding  to  one  of 
the  perturbed  Williams’  solutions.1  The  pointwise  displacements  are  then 
computed  using  the  ordinary  BIE  method.  As  in  Section  3.3.1  the  results 
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This  way  the  term  going  like  p m Ln  (p)  will  satisfy  homogeneous 
boundary  condition  along  <f>  = it. 
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where  the  «mj's  are  expansion  coefficients,  to  be  determined  and  the  expan- 
sion functions  are  the  perturbed  Williams'  solutions  which  consist  of  a 
series  of  the  Williams'  solution  forms  as  well  as  the  Williams'  logarithmic 
functions.  The  integers  M and  N are  defined  as  in  Section  3.3.  What 
results  will  be  called  the  m-th  "modified  perturbed  solution"  corresponding 
to  the  m-th  eigenfunction  of  the  infinitesimal  crack. 


Ho  insure  the  equilibrium  of  the  applied  tractions,  a solution  corres- 
ponding to  a concentrated  force  passing  through  the  center  of  the 
inner  piece  and  parallel  to  the  crack  (a  Williams'  logarithmic  solution 
with  \ = 0)  is  added  to  the  perturbed  solution. 
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The  process  is  carried  out  for  the  first  n perturbed  Williams' 
solutions.  As  mentioned  in  Section  3.3.2,  these  sets  of  solutions  tie 
the  inner  piece  to  the  annular  region  of  the  large  body. 


4.4  Evaluation  of  the  Unknown  Coefficients  of  the  Modified  Perturbed 


Solutions 


Consider  an  annular  region  of  the  body  with  the  inner  radius  fixed 
at  Rq  and  the  outer  radius  at  R (see  Fig.  10)  with  the  two  extra  regions 
ABCD  and  A'B'C'D'  added  to  it  to  make  a full  360°  (slit)  annulus.  This 
is  done  for  computational  ease  as  will  be  clear  shortly.  The  annulus  is 
loaded  along  the  boundary  with  the  traction  and  displacements  of  the  m-th 
unit  field.  The  displacement  field  for  the  inner  boundary  of  the  added 
regions  BD  and  B ' D*  will  be  defined  by  any  smooth  extrapolation.^  Along 
the  outer  boundary  L2  and  the  edges  AB  and  A'B ' , the  annulus  is  loaded 
with  the  tractions  and  displacements  of  the  modified  perturbed  solutions 
of  the  unit  field  as  given  by  Eqs.  (4.9). 

The  k-th  Williams'  solution  corresponding  to  an  infinitesimal  crack  is 
chosen  as  an  auxiliary.  The  work  integral  is  now  written  for  the  extended 
annular  region  as  follows. 


(4.10) 


r„k  ■ J „ “mj  tVdk  - V*  * °> 


♦ UJk4  = 2»)],  M > k > -N 


lThe  traction  field  is  just  the  analytic  continuation  of  the  expressions 
given  in  Eqs.  (2.12)  and  (B.7). 
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Taking  the  limit  as  R-»Rq,  we  get 

* S 

E °mj  Vjk 


mk 


j-N 


(4.11) 


where 


* /$=2  it  _ * 

'mk  =4=o  t[upm°ppk 


U4>m  ^ p k 


- ★ * . 

- [o  u . + o u..  ]}  R d$  , 
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(4.12) 


* f*=2V  r * , 

V *to  ,[dpJ  Vk ♦ D*j  >,#  3 


' [s„pj  V ♦ sp*j  V 11  V* 


(4.13) 


w 

The  can  be  evaluated  using  the  orthogonality  conditions  of 
Appendix  0. 1 


The  approximation  of  the  integral  expressions  (4.10)  for  all  k give 
us  M + N-6  simultaneous  linear  equations  with  the  a's  as  unknowns.  This 
strongly  diagonal  system  of  equations  Is  then  solved  for  the  a's.  This 
process  is  repeated  to  obtain  the  series  representation  of  each  unit  load. 
4.5.  Over-all  Solution  for  the  Loaded  Body  Having  a Wide  Crack 

Once  the  unknown  coefficients  defining  the  modified  perturbed 
solutions  are  determined,  the  general  solution  form  for  the  annulus 


The  results  of  Appendix  D are  applicable  only  because  adding  the  two 
pieces  ABCD  and  A'B'C'D'  (see  Fig.  10)  made  the  limits  of  integration 
the  same  for  the  inner  and  outer  boundaries  of  the  annulus.  This  was, 
of  course,  the  reason  for  the  extension  of  the  annulus  as  described 
earlier  in  this  section. 
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can  be  represented  by  Eqs.  (3.8)  with  the  K's  as  unknown  generalized 
displacement. 

The  solution  for  the  large  body  formed  by  the  removal  of  the  inner 
piece  is  similar  to  the  solution  for  the  case  of  the  notch  with  a fillet 
of  small  radius,  given  in  Section  3.5.  However,  the  method  will  be  modi- 
fied by  splitting  the  body  as  described  In  Appendix  A for  the  body  with  an 
infinitesimal  crack.  As  explained  in  that  Appendix,  the  straightforward 
ordinary  BIE  method  will  fall  for  the  infinitesimal  crack.  This  implies 
that  the  equations  will  be  poorly  conditioned  for  the  problem  of  the  so- 
called  wide  crack  in  a straightforward  computation  without  splitting  the 
body. 

In  order  to  simplify  the  computation  of  the  integrals  along  the 
boundary  of  the  annulus  (ABCDE  in  Fig.  11)  an  equivalent  path  of  integra- 
tion is  chosen1  (EFA  in  the  same  figure).  This  path  lies  inside  the 
corresponding  infinitesimal  crack  region  in  which  the  modified  perturbed 
solutions  are  continuous  (except  possibly  at  point  E,  the  tip  of  the 
fictitious  crack).  The  path  also  takes  into  account  the  possible  jump  of 
the  solutions  across  the  half  line  <j>  = w.  Since  some  of  the  solutions 
are  singular  at  point  E,  a piece  of  small  radius  e centered  at  E is  taken 
out  from  the  body  before  application  of  Betti's  law  to  the  large  body. 


u 

D 

Compare  Section  3.5. 


i: 
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5.  NUMERICAL  TREATMENT  OF  BOUNDARY  VALUE  PROBLEMS 
AND  SOLUTION  OF  SOME  SAMPLE  PROBLEMS 

5 . 1 Numerical  Treatment  of  Boundary  Value  Problems 

5.1.1  Preliminary  Remarks 

The  new  boundary  integral  equation  method  developed  in  Chapters 
3 and  4 reduces  the  solution  of  the  problems  under  consideration  to 
solutions  of  sets  of  singular  integral  equations.  These  equations  must 
then  be  approximated  and  the  resulting  linear  algebraic  equations  solved 
simultaneously  for  the  unknown  boundary  values. 

In  order  to  approximate  the  integral  equations,  the  closed  boundary  L is 
divided  into  N^  segments  by  N^  nodal  points.  In  the  intervals  with 
regular  boundary  points  the  tractions  and  displacements  are  specified 
at  the  nodal  points.  In  the  intervals  with  an  irregular  boundary  point, 
the  proper  Williams'  solutions  (Eqs.  2.15)  with  unknown  generalized 
displacements  represent  the  boundary  values  from  the  irregular  point  up 
to  the  so-called  transition  points.  Similarly,  in  the  fillet  and  the 
crack  problems  discussed  in  Chapter  3 and  4 the  boundary  data  on  the 
annulus  are  characterized  by  the  generalized  displacements  of  the  extended 
Williams'  solutions  (Eq.  3.7). 

The  general  form  of  the  integral  equations  for  regular  boundary 
points  is  given  by  Eq.  (2.5).  The  kernels  of  this  equation  are  the  field 
values  of  the  Kelvin  concentrated  force  solutions.  A second  type  of  integral 
equation,  given  by  Eq.  (2.17),  is  used  for  the  region  where  the  boundary  values 
are  represented  by  either  the  proper  or  the  extended  Williams'  solutions. 

The  kernels  of  this  equation  are  the  field  values  of  the  improper  Williams' 
solutions. 


IM liAriM MiH 
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The  general  form  of  the  Integrals  Involved  can  be  written  as 


I =J  G(s)ds 


G(s)  = f(s)  g(sn,s) 


(5.1) 


where  f(s)  and  g(Sg,  s)  represent  the  boundary  values  and  the  kernel 
respectively.  The  arc  length  s is  measured  from  the  singular  point  of 
the  kernel  with  s = sQ  = 0. 

The  evaluation  procedure  of  the  integrals  along  an  interval  with 
bounded  kernel  will  be  given  in  Section  5.1.2.  Expressions  for  singular 
integrals  involving  Kelvin  and  Williams'  auxiliaries  are  given  in  Sections 

5. 1.2.2  and  5. 1.2.3.  Finally  the  estimation  of  errors  is  discussed  in 
Section  5.1.3. 

5.1.2  Approximation  of  the  Integrals 
5. 1.2.1  Non-Singular  Intervals 

As  mentioned  in  the  previous  section,  the  boundary  is  divided  into 
segments  with  nodal  points.  The  length  of  the  j-th  interval  is 
denoted  by 


hj  = Sj+1  ' Sj 


(5.2) 


The  arc  length  s^  is  measured  from  the  m-th  point  (point  P in  Fig.  12), 
at  which  the  auxiliary  solution  is  unbounded.  Let  us  evaluate  the  Eq  (5.1) 
in  the  interval  [sj_i »sj+i L 


G(s )ds 


5 * tVrVil 


(5.3) 
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Since  G(s)  is  analytic  in  this  interval,  it  can  be  expanded  in  a 
Taylor  series  about  point  J 


g • gj  * $2  * 0<s3) 


(5.4) 


Now  substituting  for  (-—).  and  (^— y).  their  central  differences  and 

1,5  ^ 3s^  ^ 

performing  the  integration  we  get 

Ij  3 { [2  + a - a2]( fg) j .J  + + 3 ♦ 3a  + aZ](fg)j 


♦ [~  + 1 ♦ 2a](fg)j+i } + 0(hj  .j) 


(5.5) 


where  the  lower  case  subscripts  indicate  the  point  at  which  the  functions 

hi 

are  evaluated.  Here  a = h* — and  e = 4 for  h . f h . , , & = 5 for  h . . = h.. 

n j _ ] J J * I J - 1 J 

The  equation  (5.5)  reduces  to  Simpson's  rule  for  h.  . = h.. 

J"  ' J 

For  approximation  of  Eq.  (5.1)  in  the  interval  adjacent  to  the  singular 
point  where  there  is  a rapid  change  of  the  kernel,  a method  of  subdivision 
is  used  by  first  representing  the  density  function  as 


f(s)  * a + bs  + cs‘ 


(5.6) 


and  then  using  Simpson's  rule  for  pairs  of  adjacent  sub-intervals. 

In  the  intervals  where  the  density  function  is  represented  by  either 
the  proper  or  the  extended  Williams'  solutions,  the  kernels  are  approximated 
by  a p-th  order  polynomial.  Integration  of  Eq.  (5.1)  can  then  be  performed 
directly. 


A more  detailed  discussion  of  the  approximation  of  the  boundary  curves 
and  the  integrals  involving  Kelvin  concentrated  forces  was  given  by  Forbes[6], 
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5. 1.2.2  Intervals  Containing  the  Kelvin  Singularity 

The  value  of  the  integrals  at  the  singularity  Involving  the  Kelvin 
field  was  given  in  Section  2.2.2.  In  the  interval  containing  the  singularity, 
the  integrals  are  evaluated  by  expressing  the  boundary  values  and  the  kernels 
in  power  series  with  respect  to  the  arc  length  s measured  from  the  singular 


U 

u 

u 


point.  It  is  assumed  that  the  density  functions  are  expressible  in  the 
form  given  by  Eq.  (5.6).  The  Kelvin  kernel  can  be  expressed  as  1 

g(s)  ■ -r'+BnLn  |s|  + V a„s"  (5.7) 


♦ B0Ln  Is | *■  l ansn 
n=0 


it  is  then  possible  to  substitute  f(s)  and  g(s)  into 
.h. 


I 


nr*- 1 


m- 1 


f(s)g(s)ds 


(5.8) 


which  is  to  be  evaluated  in  the  Cauchy  principal  value  sense. 

i 

5. 1.2.3  Intervals  Containing  the  Williams1  Singularity 

The  singular  integrals  containing  the  Williams'  solution  forms  can  be 
evaluated  using  the  orthogonality  conditions  discussed  in  Appendix  C.  As 
an  example  let  us  evaluate  a singular  integral  along  the  boundary  of  a sharp 
notch  (see  Fig.  12)  with  the  n-th  Williams'  solution  form  as  the  f functions 


and  the  m-th  improper  Williams'  solution  as  the  kernel.  We  get 

* »n  * C-> 


+ X1 
n m 


HU  ■ 0)  f p 
\n  ♦ \!-i  J~c 

* HU  8 0Q)  ^ p dp 


dp 


See  Appendix  F. 
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where  the  first  term  is  the  value  of  the  integral  on  the  circular  boundary 
L£  (see  Fig.  12).  This  piece  is  taken  out  before  application  of  the  Betti's 
law  to  the  whole  body.  The  integrals  F^  and  are  given  Eqs.  (C.5)  and 
(C.8).  Here  h is  the  length  of  the  interval  between  the  irregular  point  and 
the  transition  points.  We  evaluate  1(e)  by  performing  the  integration 
and  taking  the  limit  e-*-0.  For  example,  if  both  solutions  have  zero 
traction  on  the  boundary  then  H($  = 0)  = H(<j>  = eQ)  = 0 and 


lim  1(e)  * 
e-° 


for  xn  + x*  = 0 
n m 

for  x + x*  f 0 
n m 


(5.10) 


The  other  integrals  which  arise  involving  the  Williams'  logarithmic 
solutions  can  also  be  evaluated  using  the  orthogonality  conditions  of 
Appendix  C. 


5.1.3  Errors  in  the  Approximate  Solutions 

The  errors  of  the  numerical  solutions  of  the  boundary  value  problems 
considered  in  this  study  come  from  the  following  sources. 

i.  Discretization  error  due  to  replacement  of  the  integrals  by 
finite  sums. 

ii.  Truncation  error  of  the  series  representing  the  extended  Williams' 
solutions;  and 

iii.  Round-off  errors  arising  from  evaluation  of  the  auxiliary 
fields  and  the  solution  of  the  final  linear  equations  for  the 
unknowns . 


■ ■ 
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Experience  has  shown  that  the  errors  of  the  first  kind,  the  discretiza- 
tion errors,  are  larger  than  the  other  two.  No  attempt  will  be  made  to  give 
quantitative  estimate  of  the  errors  other  than  comparing  the  result  with 
exact  solution  when  one  is  available.  In  fact,  in  all  the  problems  solved 
the  final  equations  are  strongly  diagonal  and  the  solutions  are  very  stable. 
The  solutions  do  not  change  appreciably  by  changing  either  the  quadrature 
formula  used  or  the  number  of  the  nodal  points.  The  discussion  of  the 
results  and,  in  particular,  of  these  points  will  be  given  in  Section  5.3. 

i 

5.1.4  Numerical  Implementation  ot  the  Proposed  Method  of  the  Boundary 
Value  Problems  with  Tractions  Prescribed 

We  have  seen  that  the  solution  of  a given  problem  involves  two  types 
of  integral  equations.  One  of  these  types  is  written  for  the  regular 
boundary  points  and  the  other  type  for  the  part  of  the  boundary  where  the 
solution  is  represented  either  by  the  proper  or  the  extended  Williams' 
solutions.  The  integral  equations  are  then  approximated  by  a set  of  linear 
algebraic  equations.  The  system  can  be  presented  in  matrix  form 

[A]  { u>  = {B}  (5.11) 

where  [A]  is  the  coefficient  matrix.  The  entries  of  the  column  vector  u 
are  all  the  unknowns  of  the  problem,  the  generalized  displacements,  the 
unknown  boundary  displacements  at  regular  points  and  any  unknown  tractions 
along  the  cut  (in  the  case  of  the  crack  problem).  The  vector  B is  known. 


■ 
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Before  solving  Eq.  (5.11)  for  the  unknowns,  the  rigid  body  displacements 
must  be  specified,  otherwise  the  matrix  [A]  is  singular.  This  is  done  by 
setting  three  unknown  displacements  equal  to  zero  and  eliminating  their 
corresponding  equations  from  the  Eq.  (5.11).  These  constraints  must  be 
such  that  they  prevent  rigid  body  motions  without  inhibiting  any  strains 
due  to  loading. 

Once  the  unknowns  are  found,  the  displacements  at  any  point  of  the 
body  can  be  determined  using  Eq.  (2.5).  The  determination  of  stresses 
is  discussed  in  Appendix  F.  The  field  inside  the  annular  region  can  be 
found  using  Eq.  (3.7).  The  displacements  and  stresses  at  any  point  inside 
the  inner  piece  can  be  determined  either  by  applying  the  BIE  method  to  the 
body  with  known  boundary  values  or  by  superposing  the  previously  found 
unit  fields  in  case  of  boundary  points. 

5.2  Solution  of  Some  Sample  Problems 

5.2.1  General  Remarks 

Results  of  approximate  solution  of  several  problems  using  the  integral 
equation  method  developed  in  this  study  are  given  in  the  next  sections  of 
this  Chapter.  In  all  problems,  Poisson's  ratio  taken  to  be  v = 0.25. 

Before  calculations  can  be  carried  out,  a scale  must  be  chosen  for 
the  eigenfunctions  and  other  extended  Williams'  solutions.  In  the  problems 
under  consideration,  the  magnitude  of  the  vector  displacement  for  each 
extended  Williams'  solutions  is  set  equal  to  rQ  at  e = 0 and  p = rg 
where  rQ  is  shown  in  each  figure  of  the  respective  problem. 
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The  results  presented  In  the  following  sections  Include  the  displacements 
at  certain  boundary  points  and  stresses  at  several  Interior  points.  Some  of 
the  generalized  displacements  are  also  tabulated.  All  the  values  are  given 
in  terms  of  the  applied  loadog  and  rg. 

Except  for  the  second  part  of  the  first  problem,  no  exact  analytical 
solution  exist  for  the  problems  solved  In  this  study.  However,  the  results 
of  the  second  problem  are  compared  with  Barone's  [18]  solution.  For  the 
other  problems  the  result  are  checked  by  comparing  the  solution  with  a 
sharp  notch  or  Infinitesimal  crack.  General  discussion  on  the  results 
will  be  given  in  Section  5.3. 

5.2.2  A Cracked  Disk  Under  Constant  Radial  Tension  (Problem  No.  1) 

As  discussed  in  Section  2.3  the  main  object  of  this  example  is  to 
demonstrate  how  an  arbitrary  self-equilibrated  load  on  the  disk  can  be 
expanded  using  proper  W1 11 iams'  solutions.  The  disk  as  shown  in  Fig.  16 
is  loaded  with  radial  traction  of  constant  magnitude  T = const.  Equations 
(2.15)  are  used  for  the  expansion  of  the  load.  The  improper  Williams' 
solutions  corresponding  to  the  infinitesimal  crack  (Eqs.  2.16)  are  taken 
as  auxiliaries.  Applying  the  reciprocal  theorm  to  the  body  in  which  the 
singular  point  P is  excluded  by  a disk  of  small  radius  e and  then  taking 
the  1 imit  e+0,  we  get 

M xm+xd+m2*  * * 

m l 4 Vo  Jj4,0  tV*>  PJ<*>  * 

m / j 

♦oj j*  ” [PjU)  A*U)  + RjU)  B*U)]d*  (5.11) 

X%l/*=2n  * * 
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where  the  a's  are  the  unknown  expansion  coefficients,  ^is  the  radius  of  the 
of  the  disk,  and  M is  total  number  of  expansion  function  including  the  three 
rigid  body  motions  (supressed  in  this  problem). 

The  left  hand  sides  of  the  integral  equations  were  computed  directly 
and  the  right  hand  sides  were  approximated  by  Simpson's  rule.  Table  1 
gives  the  expansion  coefficients  for  three  different  numbers  of  eigen- 
functions used  as  terms  in  the  expansion.  It  can  be  seen  that  the  values  of 
the  coefficients  stabilize  as  the  number  of  functions  increases.  This  is 
especially  true  for  the  first  coefficient  which  is  proportional  to  the  stress 
intensity  factor.  Table  2 gives  the  computed  value  of  tractions  on  the 
boundary.  Convergence  to  the  exact  answers  depends  on  the  character  of  the 
applied  tractions.  For  example,  this  same  problem  was  solved  with  loading 
not  uniform  but  corresponding  to  the  Griffith  problem  (Fig.  17).  The 
results  are  given  in  Tables  1 and  2 along  with  those  of  constant  radial 
loading.  It  can  be  seen  that  convergence  is  much  better  than  for  the  first 
loading  case.  This  is  expected  because  the  first  few  Williams'  solutions 
do  represent  the  field  in  the  vicinity  of  the  crack.  This  will  be  the 
case  in  all  actual  calculations  in  real  stress  concentration  problems  pro- 
vided that  the  circle  is  reasonably  small  with  the  body  diameter  and 
greater  than  the  fillet  radius. 

As  was  mentioned  the  first  coefficient  of  the  expansions  in  Table  1 
is  proportional  to  the  stress  intensity  factor.  For  example,  for  the 
Griffith  load  [33]  and  M = 17  this  value  is  to  within  0.04  percent  of  the 


exact  value. 
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5.2.3  Uniformly  Loaded  Plate  with  a Straight  Edge  Crack  (Problem  No.  2) 

This  problem  (Figs.  18  and  19)  Is  solved  by  "splitting"  the  body  as 
described  in  Appendix  A.  In  this  method,  the  crack  is  extended  by  an 
imaginary  cut  (the  dashed  line  AB  in  Fig.  18);  then  the  Barone  and  Robinson 
solution  technique  is  applied  to  both  pieces  of  the  body  simultaneously. 

The  generalized  displacements  for  the  irregular  point  A (point  1 in 
Fig.  19)  are  given  in  Table  3.  The  displacements  for  some  boundary  points 
are  given  in  Table  4,  along  with  the  values  from  Barone's  [18]  solution. 

As  can  be  seen  they  are  in  good  agreement.  The  stress  distribution  along 
the  cut  AB  is  given  in  Fig.  20.  Although  the  stresses  resulting  from  the 
Barone  solution  are  not  given  here,  they  also  match. 

It  should  be  mentioned  that  the  introduction  of  the  Williams'  solution 
at  point  25  at  the  end  of  the  cut  (Fig.  19)  is  for  computational  ease  and 
reduction  of  the  number  of  unknowns;  this  point  could  have  been  treated 
as  a regular  point. 

5.2.4  Uniformly  Loaded  Plate  with  a Curved  Edge  Crack  (Problem  No.  3) 

This  problem  (see  Figs.  21,  22  and  23)  is  solved  to  show  the  versatility 
of  the  "splitting"  method  described  in  Appendix  A.  The  solution  procedure 
is  essentially  the  same  as  that  for  the  straight  edge  crack  (Problem  No.  2). 
The  generalized  displacements  for  the  tip  of  the  crack  is  given  in  Table  3 
along  with  those  of  problem  2 (straight  edge  crack).  Here,  the  generalized 
displacements  of  mode  5 and  6 can  be  used  for  the  determination  of  the 
stress  intensity  factor  for  normal  and  shearing  mode. 

The  displacements  for  some  boundary  points  are  given  in  Table  5.  The 


distribution  of  the  normal  and  radial  stresses  on  the  cut  are  given  in  Fig. 
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24.  Curiously  enough,  the  normal  stress  distribution  on  the  particular  cut 
chosen  is  the  same  as  for  the  straight  edge  crack1  and  is  not  a function 
of  the  shape  and  orientation  of  the  cut  so  long  as  the  crack  remains  inside 
the  left  hand  side  edge  of  the  body. 


5.2.5  Uniformly  Loaded  Body  with  a Notch  Having  a Fillet  of  Small  Radius 
(Problem  No.  4) 

This  problem  as  shown  in  Figs.  25,  26  and  27  is  solved  using  the  method 

\ 

described  in  Chapter  3.  The  field  in  the  annular  region  is  represented 
2 

by  four  real  and  four  complex  modified  Williams'  solutions.  Since  each 

\ 

complex  solution  gives  two  independent  real  solutions  [Ref.  18],  we  have 
twelve  generalized  displacements  as  unknowns  for  the  annular  region. 

\ 

v Table  6 gives  the  expansion  coefficients  for  the  first  three  modified 
Williams'  solutions  using  nine  proper  and  three  improper  eigenfunctions. 

As  described  in  Section  3.4  these  solutions  are  derived  from  the  expansion 

\ 

of  the  unit  fields.  Table  7 gives  one  component  of  displatements  and 
tractions  of  the  first  two  unit  fields  along  with  their  corresponding 
values  from  the  modified  Williams'  solutions.  It  can  be  seen  from  the  table 
they  are  remarkably  close.  Table  8 gives  the  generalized  displacements  of 
the  modified  Williams'  solutions  for  the  annular  region  (see  Fig.  4). 

The  displacements  for  some  boundary  points  are  given  in  Table  9.  The 
same  table  gives  the  results  from  the  solution  of  the  fillet  problem  using 


From  the  geomentry  and  the  loading  of  the  two  problems,  the  integral  as 
well  as  the  first  moment  of  the  normal  stress  along  the  cut  should  match 
for  obvious  statical  reasons,  but  it  is  interesting  that  their  distribu- 
tions also  match. 


2 

Including  two  rigid  body  translations  and  one  rotation. 


the  ordinary  BIE  method  with  50  boundary  points.  The  stresses  and  displace- 
ments for  certain  points  of  the  inner  piece  are  given  in  Table  10.  Finally, 
the  stress  distribution  along  the  axis  of  symnetry  is  given  in  Fig.  28. 

5.2.6  Uniformly  Loaded  Body  with  a Wide  Edge  Crack  (Problem  No.  5) 

A uniformly  loaded  body  with  a parallel-sided  crack  of  finite  but 
small  width  (wide  crack)  as  shown  in  Figs.  29,  30  and  31  is  solved  by  the 
method  described  in  Chapter  4. 

As  discussed  in  Section  4.2,  the  first  step  of  the  solution  is  deter- 
mination of  the  perturbed  Williams'  solutions.  Tables  11  and  12  give  the 
stresses  and  displacement  of  the  perturbed  and  unperturbed  Williams' 
solution  for  two  points  at  the  boundary  of  the  inner  piece.  As  can  be 
seen  from  the  tables,  the  perturbation  alters  the  magnitude  of  the  displace- 
ment field  of  the  unperturbed  solution  by  less  than  5 percent  (for  the 
first  solution).  The  effect  of  perturbation  on  the  stress  field  Is  less 
pronounced  for  points  inside  the  body  and  away  from  the  boundary.  Table  13 
gives  expansion  coefficients  for  3 modified  perturbed  solutions  and 
Table  14  gives  one  component  of  traction  and  displacement  to  two  unit 
fields  and  their  matching  values  from  the  expansion  (the  modified  perturbed 
solutions).  The  results  are  not  as  good  as  in  the  corresponding  expansion 
for  the  fillet  problem  (Section  3.3),  especially  for  the  traction  in 
x-dlrection. 

The  generalized  displacements  for  the  annulus  are  given  in  Table  15 
and  displacements  for  a selection  of  boundary  points  are  shown  In  Table  16 
and  finally  the  stress  distribution  along  the  cut  is  given  in  Fig.  32. 
Comparison  of  the  stress  distribution  with  that  of  the  infinitesimal 
crack  (Fig.  20)  shows  that  except  for  the  region  in  the  vicinity  of  the 
crack  they  almost  coincide,  which  is  expected. 
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5 . 3 Discussion  of  the  Results 

The  solutions  presented  in  the  previous  sections  were  found  to  be 
stable  with  respect  to  the  number  of  boundary  points  and  the  number  of 
extended  Williams'  solutions  used  in  the  expansion  expressions.  This  is 
also  the  general  conclusion  reached  by  Barone  [18]  for  the  infinitesimal 
crack  and  sharp  notch  problems.  This  stability  reflects  the  fact  that 
the  proposed  integral  equation  approach  leads  to  a strongly  diagonal 
system  of  linear  algebraic  equations.  The  accuracy  of  the  integral 
approximation  is  also  an  important  factor.  Certain  numerical  refinements 
such  as  subdividing  the  boundary  intervals  and  including  the  boundary 
curvature  in  the  evaluation  of  the  singular  integrals  played  a particularly 
important  part  in  achieving  accurate  results. 

In  the  case  of  the  wide  crack  problem,  the  expansion  of  the  first 
unit  field  is  not  as  accurate  as  was  expected.  It  seems  that  for  more 
accurate  expansion,  several  more  perturbed  Williams'  solutions  would  be 
needed.  However,  each  additional  Williams'  solution  involves  traction 
that  vary  rapidly  with  <j>  on  the  circular  boundary.  Stresses  at  the 
curved  tip  of  the  crack  have  been  found  applying  such  variation  to  the 
boundary  the  stress  which  arise  are  very  small  (less  than  0.25  percent 
of  the  stress  applied  to  the  boundary).  In  addition,  calculation  was 
done  applying  the  difference  between  the  applied  traction  of  the  unit 
field  and  the  traction  found  from  the  truncated  series.  Again  no  significant 
stress  is  found  at  the  roots  of  the  crack.  It  is  therefore  unnecessary  to 
refine  the  expansion  for  the  annulus  because  the  stress  pattern  does  not 
change  where  it  is  of  interest. 

I I 


m 

a I 


51 


6.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDIES 


6.1  S unwary 

A boundary  Integral  equation  technique  has  been  described  for  the 
analysis  of  homogeneous  elastic  bodies  for  which  small  displacement  theory 
Is  assumed  to  hold.  Major  emphasis  Is  placed  on  treatment  of  bodies  with 
high  stress  concentrations.  Examples  of  such  problems  are  bodies  with 
notches  having  a fillet  of  small  radius  and  bodies  with  parallel-sided 
cracks  of  small  but  finite  width.  The  formulation  of  the  solution  technique 
Is  based  on  the  reciprocal  theorem.  The  boundary  Integral  equation  technique 
uses  only  the  boundary  values  not  the  Interior  field  for  the  solution  of  a 
given  problem. 

The  method  of  solution  for  the  notch  with  a fillet  of  small  radius 
is  developed  from  the  method  Introduced  by  Barone  and  Robinson  [18],  for 
the  solution  of  a sharp  notch.  At  moderate  distance  from  the  fillet,  the 
intermediate  field  In  the  body  Is  expressible  In  terms  of  solutions  of 
Williams'  type  [16].  Unlike  the  case  of  the  sharp  notch,  decaying  as  well 
as  the  growing  solutions  are  present.  An  Inner  region  containing  the  fillet 
Is  separated  out  and  analyzed  under  load  systems  corresponding  to  the  Inter- 
mediate fields  of  Williams'  type.  It  Is  then  possible  to  match  the  displace- 
ments and  stresses  on  the  boundary  between  the  Inner  and  Intermediate  bodies. 
It  Is  In  this  way  that  the  details  of  the  boundary  configurations  of  the 
fillet  are  reflected  In  a set  of  generalized  coordinates  which  characterize 


l 
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the  Intermediate  field.  The  remainder  of  the  problem,  for  the  intermediate 
and  the  distant  field,  is  solved  in  essentially  the  same  way  as  Barone  and 
Robinson  [18]  handled  the  sharp  notch. 

The  solution  of  the  wide  crack  problem  is  similar  except  that  it  is 
necessary  to  alter  the  Williams'  solution  applicable  to  the  infinitesimal 
crack.  Perturbed  Williams'  solution  are  developed  for  this  purpose. 

Finally,  a modified  version  of  the  Barone  and  Robinson  solution  technique 
is  used  to  avoid  certain  numerical  difficulties.  The  modification, 
incidentally,  is  very  useful  for  the  infinitesimal  crack  problem  as  well. 

6.2  Conclusions 

The  general  conclusion  of  this  study  is  that  the  method  developed 
gives  accurate  and  numerically  stable  results  with  rather  small  computational 
effort.  The  solution  of  the  problems  is  fairly  insensitive  to  the  number 
of  extended  Williams'  solutions  used  for  the  region  with  high  concentration 
of  stresses.  The  solution  also  Is  not  affected  very  much  by  the  number  of 
boundary  points  used.  Instead,  It  depends  more  on  the  method  of  approxi- 
mation used  for  the  evaluation  of  the  integrals.  In  particular,  a larger 
interval  length  can  be  used  provided  there  is  finer  subdivision  for  the 
evaluation  of  the  Integrals.  However,  no  quantitative  analysis  has  been 
carried  out  to  determine  the  optimum  interval  length  and  minimum  computing 
time. 

The  methods  of  this  stucty  can  be  used  very  effectively  for  determination 
of  the  stress  intensity  or  stress  concentration  factor.  For  example,  the 
stress  Intensity  factor  for  the  Griffith  problem  [33]  was  computed  to 
within  0.3  per  cent  of  the  exact  solution  using  the  splitting  method  of 
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Appendix  A.  In  this  computation  there  were  30  boundary  points  with  51 

unknowns  (taking  advantage  of  the  symmetry)- 

/ 

6.3  Recomnendations  for  Further  Studies 

One  important  extension  of  the  methods  developed  in  this  study  would 
be  an  analysis  in  which  the  material  is  treated  as  elastic-plastic  rather 
than  elastic.  The  application  of  load  to  bodies  with  sharp  notches  will 
often  cause  formation  of  a plastic  region  limited  to  the  vicinity  of  the 
singular  point.  This  is  also  true  for  bodies  with  high  stress  concentration 
such  as  those  considered  in  this  study.  A method  similar  to  that  used  for 
the  solution  of  the  fillet  problem  could  be  used  for  this  case.  The  body 
would  be  divided  into  three  regions,  an  inner  piece  containing  the  fillet 
and  the  entire  plastic  region,  an  annular  region  for  which  a modified 
Williams’  solution  with  non-linear  coefficients  would  be  developed,  and  the 
larger  body  which  is  tied  to  the  annulus  by  the  usual  integral  equation 
technique. 

Another  problem  that  deserves  consideration  is  the  solution  of 
piecewise  non-homogeneous  bodies  with  a crack.  This  problem,  as  well  as 
the  problem  of  non-homogeneous  bodies  with  a V-notch  may  be  solved  with 
no  apparent  difficulty  using  the  "splitting"  technique  discussed  in 
Appendix  A.  The  stress  and  displacement  field  in  the  vicinity  of  the  crack 
tip  that  would  be  needed  for  the  solution  of  the  problem  by  this  method,  is 
given  by  O.K.  Aksentian  (24). 

A third  problem  of  interest  is  the  improvement  of  the  perturbed 
Williams'  solutions  developed  in  Chapter  4 for  the  solution  of  the  wide 
crack  problem.  The  perturbation  procedure  developed  seems  to  have  an 
asymptotic  rather  than  a convergent  character.  If  the  opening  angle  of 
the  inner  piece  (angle  DFB  in  Fig.  11)  is  sizable,  the  accuracy  by  this 
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perturbation  method  may  be  limited  to  the  degree  where  it  is  unsatisfactory. 
The  development  of  a convergent  procedure  would  be  very  desirable.  Several 
unsuccessful  attempts  were  made  to  find  a simple  procedure  in  the  course 
of  analysis. 

Another  extension  would  be  to  develop  a numerical  method  for 
evaluating  the  effect  of  the  plate  thickness  on  the  singular  stress 
distribution.  Analysis  of  thick  plates  shows  that  the  stress  state  is 
essentially  three-dimensional  near  the  edge  [34-36].  However,  the  fact 
that  this  boundary  layer  disturbance  is  local  suggests  the  possibility 
of  applying  a local  correction  to  the  two-dimensional  solution. 
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TABLE  4 

SELECTED  BOUNDARY  DISPLACEMENTS 
FOR  PROBLEM  2 


Boundary 

Points 

Splitting  Method 
(Appendix  A) 

Barone's 

(Ref. 

i 

Solution 

is)  ; 

: Ux/a] 

Uv/a 

U /a 

X 

VB 

1 

-47.67 

| 

0.00 

-47.62 

0.00 

2 

-48.22 

-21.11 

-48.18 

-21.10 

3 

-48.65 

-31.86 

-48.60 

-31 .87 

4 

-48.85 

-41.63 

-48.80 

-41.60 

5 

-48.89 

-50.85 

-48.83 

-50.81 

6 

-48.83 

-59.70 

-48.76 

-59.61 

7 

-48.83 

-68.42 

-48.76 

-68.30 

8 

-40.14 

-68.42 

-40.07 

-68.30 

10 

! -23.44 

-68.59 

-23.43 

-68.48 

18 

-10.97 

-11.07 

-10.97 

-11.06 

21 

-16.20 

3.65 

-16.22 

3.60 

23 

-28.10 

5.10 

-26. 762 

5.13 

25 

-35.15 

0.00 

-35.06 

0.00 

1 °l 

a = j 

- rQ  where  oQ  is 

the  applied  load,  E is  the  modulus 

of  elasticity  and 

rQ  is  given  in  Fig.  19. 

2 

There  may  be  a typographical  error  in  Ref  18. 


The  first  three  modes  corresponding  to  rigi.d  body  motions  are  not  included  in  the  expansion. 
xn  is  corresponding  eigenvalue  of  the  applied  traction  used  to  develop  the  unit  field. 


TABLE  8 

GENERALIZED  DISPLACEMENTS  FOR  PROBLEM  4 


u 

u 

D 

0 

D 

0 

D 


Mode 

No. 

Eignevalue 

x1 

(General ized 

2 

Displacements) /a 

1 

0.00 

0.000 

2 

0.00 

6.719 

3 

1.00 

0.000 

4 

0. 54448 

6.140 

5 

1 . 6293+i 0.2313 

4.767 

6 

-0.274 

7 

2 . 97 18+i 0.3739 

-2.169 

8 

-1.169 

9 

4. 3104+i0.4555 

-0.338 

10 

-0.029 

11 

5. 64712+iO. 5137 

-0.535 

12 

0.326 

* X corresponds  to  the  applied  traction  of  the  unit  field. 

2 °0 

a = where  aQ  is  the  applied  load  and  E is  the  modulus  of  elasticity. 
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TABLE  9 

SELECTED  BOUNDARY  DISPLACEMENTS 
FOR  PROBLEM  4 


Point 

No.1 2 

(Displacements)/^ 

Method  of  Chapter  3 

Ordinary  BIE  Method 

Ux 

_ ", 

Ux 

c: 

1 

0.000 

6.331 

0.000 

6.330 

2 

-3.015 

5.664 

-3.018 

5.635 

3 

-4.658 

4.897 

-4.656 

4.887 

4 

-5.827 

4.153 

-5.829 

4.149 

5 

-6.578 

3.516 

-6.581 

3.510 

9 

-8.996 

-0.699 

-8.999 

-0.700 

11 

-6.533 

0.000 

-6.533 

0.000 

13 

-4.214 

0.778 

-4.214 

0.776 

17 

-0.593 

4.726 

-0.592 

4.718 

19 

0.000 

6.025 

0.000 

6.017 

1 °0 

a = Tg  where  oQ  is  the  applied  load,  E is  Young  modulus  of 
elasticity  and  rQ  is  shown  in  Fig.  27. 

2See  Fig.  27. 
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TABLE  10 

STRESSES  AND  DISPLACEMENTS  AT  SELECTED  POINTS 
OF  THE  INNER  PIECE  (Problem  4) 


1. 

Values  from: 

2. 

Method  of  Chapter  3 

Ordinary  BIE  Method 

Stresses 

Displacements 

KOI  n t 

No. 

warn ibb 

°xy/o0 

U /a3 

X 

MBM 

i 

1.363 

1.363 

1.363 

-1.962 

6.002 

1.372 

1.374 

1.363 

-1.960 

5.998 

3.464 

0.520 

-0.764 

-1.123 

5.823 

10 

3.459 

0.526 

-0.772 

-1.120 

5.820 

2.842 

1.632 

0.000 

0.000 

6.181 

28 

2.842 

1.628 

0.000 

0.000 

6.180 

13.522 

0.002 

0.000 

0.000 

6.265 

64 

13.425 

0.144 

0.000 

0.000 

6.259 

1 

7.430 

2.219 

0.000 

0.000 

6.331 

73 

7.408 

■ , . 

2.206 

0.000 

0.000 

6.330 

See  Fig.  26. 

oq  is  the  applied  load. 

°0 

a = r.  where  E is  modulus  of  elasticity  and  r_  is  shown  in  Fig.  27. 
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TABLE  13 

COEFFICIENTS  OF  EXPANSION  FOR  SELECTED  MODIFIED 
PERTURBED  SOLUTIONS  (Problem  5) 


Eigenvalues 

,1 


Coeff.  of  Expansion 


.5 

A 

= 1.50 

X 

1.06473 

0.04392 

-0.13221 

0.01062 

0.04513 

-0.03285 

0.04673 

-0.01040 

-0.03833 

-0.00039 

-0.00715 

0.00122 

-0.00086 


-0.02285 

-0.00507 

1.04189 

-0.01259 

0.01923 

-0.00215 

-0.00964 

0.01021 

-0.00669 

0.00004 

0.00158 

0.00000 

-0.00001 


= 2.00 


-0.01272 

-0.00960 

0.02557 

1.00871 

-0.01630 

0.00092 

0.02006 

-0.02406 

0.01315 

0.00123 

0.00231 

-0.00019 

-0.00008 


X corresponds  to  eigenfunctions  of  infinitesimal  crack. 

) 

"This  is  the  concentrated  force  solution  which  is  used  to  satisfy  the 
equilibrium  in  x-direction. 


TABLE  15 

GENERALIZED  DISPLACEMENTS  FOR  PROBLEM  5 


Mode 

No. 

X1 

(Generalized  Displ.)/a2 

1 

0.00 

0.000 

2 

0.00 

-47.420 

3 

1.00 

0.000 

4 

0.50 

-18.167 

5 

1.00 

0.041 

6 

1.50 

-0.753 

7 

2.00 

0.109 

8 

2.50 

0.703 

9 

3.00 

-0.308 

10 

3.50 

0.422 

11 

4.00 

-0.276 

12 

4.50 



-0.281 

* X corresponds  to  the  tractions  used  for  the  development  of  the 
unit  field. 

2 °0 

a = -jj-  where  Oq  is  the  applied  load  and  E modulus  of  elasticity. 
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TABLE  16 

SELECTED  BOUNDARY  DISPLACEMENTS  FOR  PROBLEM  5 


Boundary 

(DisplacementsJ/a1 

Points 

U 

X 

U 

y 

1 

-47.42 

0.00 

2 

-47.97 

-21.32 

3 

-49.03 

-32.29 

4 

-49.30 

-42.28 

5 

-49.38 

-51.64 

6 

-49.29 

-60.69' 

7 

-49.29 

-69.52 

8 

-40.46 

-69.52 

10 

-23.88 

-69.70 

18 

-11.19 

-11.27 

21 

-16.52 

3.68 

23 

-28.57 

5.06 

25 

-35.69 

0.00 

a = rO  where  Og  is  the  applied  traction,  E is  the  Young 
modulus  of  elasticity  and  rn  is  shown  in  Fig.  31. 


FIG.  14.  AN  ANNULAR  SECTION  OF  A BODY  WITH  A NOTCH  FOR  DEVELOPMENT 
OF  THE  ORTHOGONALITY  CONDITIONS 


B 


FIG.  15.  DETAIL  OF  THE  ANNULAR  REGION  FOR  DEVELOPMENT  OF 
ORTHOGONALITY  CONDITIONS 


r* 

O Irregular  Point 
o Transition  Point 

86 

. r0  . 

6 at 

1-  2/3r0 

ro  . 

FIG.  20.  COMPUTED  STRESSES  FOR  SECTION  AB  (Problem  2) 


FIG.  24.  COMPUTED  STRESSES  ALONG  THE  CUT  AB  (Problem  3) 


OF  SMALL 


G.  26.  DETAIL  OF 
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APPENDIX  A 

SOLUTION  OF  A CRACK  PROBLEM  BY  SPLITTING  THE  BODY 


A.l  Introduction 


Some  difficulties  arise  when  solving  a crack  problem  using  the  straight- 
forward Barone  and  Robinson  method.  The  source  of  the  problem  can  be 
explained  as  follows:  Application  of  the  Kelvin  concentrated  force 
solution  as  an  auxiliary  at  points  along  the  crack  boundary  (beyond 
transition  points),  leads  to  an  indeterminate  integral  equation,  because 
this  force  picks  out  only  two  independent  combination  of  four  displacements, 
namely,  the  sum  of  them  at  the  top  and  bottom  of  the  crack.  In  the  Barone 
and  Robinson  solution  technique  [18],  in  order  to  pick  out  the  difference 
of  the  displacements,  Westergaard's  solutions  [26]  were  used  as  an  auxiliary. 
This  solution  corresponds  to  equal  and  opposite  prying  forces  acting  at 
some  point  along  a semi-infinite  crack  in  an  infinite  body.  However, 
this  solution  is  applicable  only  for  a straight  crack.  Moreover,  the 
use  of  this  solution  represents  a very  special  technique  which  does 
not  fit  comfortably  in  the  framework  of  the  general  procedure. 

In  some  other  investigations  the  problem  has  been  resolved  by  either 
developing  a Green's  type  function  to  fit  the  entire  physical  problem  [13], 
or  replacing  the  crack  by  an  open  notch  [14].  The  first  method  is  naturally 
limited  to  very  special  types  of  crack  problems.  The  second  approach, 
which  actually  changes  the  geometry  of  the  body,  may  provide  a fairly 
accurate  value  of  the  stress  intensity  factor,  but  does  not  give  accurate 
value  of  stress  field  near  the  crack  tip. 


Li 
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A . 2 Splitting  the  Botty 

In  order  to  overcome  the  above  mentioned  difficulty,  it  is  convenient 
to  cut  the  cracked  body  by  extending  the  crack  in  some  arbitrary  fashion. 
The  Barone  and  Robinson  solution  technique  may  then  be  applied  to  each 
piece  of  the  body-  The  resultant  equations  can  be  solved  simultaneously 
along  with  constraint  equations  which  must  be  Introduced  to  insure 
compatibility  of  displacements  and  tractions  along  the  imaginary  cut. 

The  solution  to  two  sample  problems  using  the  method  are  given  in 
Section  5.2.  It  is  worth  mentioning  that  the  proper  Wi 11  lams ' solutions 
are  used  to  represent  the  stress  state  along  the  interval  extending  from 
the  irregular  point  up  to  the  so-called  transition  points.  Therefore, 
in  order  to  reflect  the  actual  nature  of  the  singularity  in  the  vicinity  of 
the  irregular  point,  the  interior  angle  from  boundary  to  boundary  (wedge 
angle)  must  be  taken  as  2n.  Even  though  each  of  the  pieces  of  the  cut 
body  has  an  angle  less  than  2 w,  the  cut  is  wholly  fictitious  and  the 
stress  field  of  interest  is  that  for  the  uncut  body. 


I p 

...  Ik 
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APPENDIX  B 

DERIVATION  OF  NEW  SOLUTIONS  FROM  THE  WILLIAMS* 
DISPLACEMENT  FIELD  (WILLIAMS*  LOGARITHMIC  SOLUTIONS) 

A new  solution  can  be  derived  by  differentiating  the  Williams’ 
displacement  field  with  respect  to  the  parameter  X,  the  exponent  of  p 
in  Eqs.  (2.10).  To  show  this,  we  rewrite  Eqs.  (2.1) 

Up  = A(*)pA 
= BU)PX 


This  displacement  field  satisfies  the  equilibrium  equations  for 
all  values  of  x,  so  does  any  linear  combination  of  these  fields  for 
different  values  of  X.  In  particular 


(B.2) 


also  satisfy  the  equilibrium  equations.  Taking  the  limit  as  x^  -*■  x^ 
we  get  a new  solution  of  the  elasticity  equations: 


(1) 

3U 

3AU) 

up 

p - 

ax 

[ 

ax 

+ A(<j>)  Ln(p )]  p 

(1) 

3U. 

= £ = 

ax 

aBU) 

% 

[ 

ax 

+ B(<f>)  Ln(p)]  p 

u 
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where 

= [cos(x-l)*  - (x-3+4v)$  sin(x-l)*]C1 

+ [sin(x-lH  + (x-3+4v)4>  cos(x-l )*]C2  (B.4) 

+ [~$  s1n( x+1  )^]Cj  + [«j>  cos(x+l  )4»]C^  , 

3 C-s1n(x-l  - ( x*3-4v)|  cos  ( x-1  U]Cj 

+ C cos ( x- 1 ) 4*  - (x+3-4v)tf  si n( x- 1 )^]C2  (B . 5 ) 

+ [-<1  cos( x+1  HK3  + [-*  sin(x*1  )^]C4 


The  corresponding  stresses  are 


(1) 

ao 

aPU) 

__H£_ 

- [ 

— 

+ P(t)Ln(p)]p 

pp 

ax 

ax 

(1) 

ao  . 

aQU) 

♦♦ 

_>i£ 

ax 

* [ 

ax 

+ Q(+)  Ln(p)]  p 

(1) 

ao  . 

aRU) 

3 

p<J> 

_J 1$L 

ax 

- [ 

ax 

+ R(<f>)  Ln(p)]  p 

where  PU),  QU)  and  R($)  are  given  in  Eq.  (2.13). 

Other  new  solutions  can  be  found  by  taking  higher  derivatives  of  the 
Williams'  displacement  fields  with  respect  to  X.  The  general  form  of  the 
n-th  solution  is  given  in  Appendix  D using  complex  fields  and  complex 
variables. 

.1 

J 


APPENDIX  C 


ORTHOGONALITY  CONDITIONS  FOR  THE  WILLIAMS'  SOLUTION  FORMS 


C.l  Orthogonality  Conditions  for  the  Extended  Wi  Haims'  Solutions 


Consider  a body  with  a sharp  notch  of  Interior  angle  0q  (see  Fig. 


14,  15).  We  remove  the  annulus  ABCDEF  from  the  body  as  shown.  Let  us 


write  Betti's  law  for  the  piece  using  two  different  Williams'  solution 


forms. 


W = W.  + W.  = 0; 
L1  l2 


(C.l) 


W.  = / [t.u.  - u*t.' 

Li  ;L  j j j j- 


1-1.2  ; J ■ 1.2 


(C.2) 


In  Equation  (C.2)  the  sumatlon  convention  is  used  for  the  lower-case 


Indices.  The  L1  is  the  curved  part  of  the  boundary  and  L?  consists  of 
the  two  straight  edges  shown  in  Fig.  15. 


Substituting  for  the  tractions  and  displacement  in  Eq.  (C.2)  from 
Eqs.  (2.15),  we  get 


Wi  = F [(p  + dp) 

L^  mnLV^ 


* * 

xm  + X X + X 

m n m n, 
- p J 


(C.  3) 


w,  * f r (x  + xjp 

mn1-  m n 


* xm  + xn-l 

m n 


(C.4) 


I *aQo 

Fmj4m0  { ty^Pn^  + Bm^ 


(C.5) 


-[A*U)  Pm(«)  + B*U)  Rn ( 4> ) ]> d<». 


n,T'  n' 


m 


l 
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xm  + x -1 

ii  f*  ^ n i 
W,  s G p dp 

mn 


G “ H (41  3 6rv)  ■ H (6  = 0) 
nn  mn'^  0'  mn  ' 


Here 


H(6)  = [Am(6)  R*U)  ♦ Bm(6)  Q*U)] 


- [y*)  1^4)  + Bn(6)  Qm(6)] 


Substituting  for  W,  and  W,  in  Eq.  (C . 1 ) we  obtain 
L1  l2 

* 

* Xm  + Xr,_1 

r( A + x ) F + G ] P n dp  = 0 
LV  m n7  mn  mnJ  ^ 


Since  Equation  (C.9)  must  be  true  for  all  values  of  p we  have 


( X + Xj  F + Gmn  = 0 
m n mn  mn 


We  consider  the  two  possible  cases: 


a)  G = 0 then  F = 0 if  Xm  + Xn  t 0 
mn  mn  m n 


(C.  6) 
(C.  7) 


(C.8) 


(C.9) 


(C.10) 


(C.ll) 


and  for  Xm  + \n  = 0 Eq.  (C.10)  does  not  give  any  information  about  Fmn> 


m n 


Conversely,  if  the  first  term  in  Eq.  (C.10)  is  zero  then  Gmn  = 0. 


b)  Then  F 


mn 


mn  x +X_ 
m n 


V for  m 


+ xn  t 0 (C.12) 


j 


i 
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which  shows  that  the  integrals  around  the  curved  boundary  can  be  evaluated 


by  finding  the  work  terms  at  the  two  end  points  and  divide  it  by  + xp. 


C . 2 Orthogonality,  Conditions  for  the  Williams'  Logarithmic  Solutions 


Just  as  in  the  previous  case,  we  write  Betti's  law  for  the  annulus 
(Fig.  15)  using  two  solutions,  where  one  of  them  is  a Williams'  logarithmic 
solution  and  the  other  one  is  an  ordinary  Williams'  solution.  We  have 


(1)  (1)  (1) 


W.  + W.  =0 
L1  L2 


(C.13) 


where  the  superscript  indicates  that  the  Williams'  logarithmic  solution 


used  is  the  first  derivative  of  the  Williams'  solution  with  respect  to 


,(D 


the  exponent  of  p.  Evaluating  W,  ' using  Eq.  (C.2),  we  get 

L1 


★ * 

\ + X X + X 

m n m n-, 
- p J 


+ Fmn  [Ln(p  + dp)  (p  + dp)  m n -Ln(p) 


* 

Am  + X„ 
m nn 

p J 


(C. 14) 


or 


(1)  * x_  + x -l 

W,  = F (X  + X ) p 01  n 

1 mn  ' m n • ► 


mn  ' m n' 


dp 


(C. 15) 


* p ♦ i*.  * 0 n<p>k v' 


dp 
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where  Fmn  Is  given  by  Eq.  (C.5)  and 


m A3  en  (1)  * (1) 

» 'A-o  ltAm  <♦>  \(*>  ♦ 

* 0)  * 0) 

- [AnU)  PmU)  + BnU)  R^UHld* 


(C. 16) 


(1)  , . 
The  W.  Is  also  evaluated  using  Eq.  (C.2): 
l2 

(!)  A+  A*-l  (1)  Xm  + V1 

WU  &m  0 dp  * Gmn  ’ 


Ln(p)dp 


(C. 1 7) 


where  G„„  is  given  by  Eq.  (C.7)  and 
mn 


Oi  0)  u> 

Gnr  ’ % (♦  ’ e0>  - Hmn  <♦  ’ 0)i 

(1)  0)  * (1)  * 

Hn,n  <♦>  ■ * Bm(»)  <)„(♦» 

* (1)  * (1) 

- C^U)  Rm(«)  * BnU)  Qm(+)] 


0)  0)  , , 

Substituting  for  W.  and  W.  In  Eq.  (C.19),  we  get 
L1  l2 


* (1)  S'K 

+ xn^  Fmn  + Pmn  + ^mn 


* + a-1 

♦ t(\n  + O L ♦ O Ln(P)}  e n dp  3 o 

*■  m n mn  mn 


(C.18) 


(C.19) 


(C. 20) 


Since  Eq.  (C.20)  must  be  true  for  all  values  of  p,  we  obtain 


i ’ M>dB 
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APPENDIX  D 

THE  WILLIAMS'  SOLUTIONS  IN  COMPLEX  FORM 

Most  calculations  using  Williams'  solutions  are  carried  out  in  the 
natural  polar  coordinates.  However,  for  the  wide  crack  problem.  It  is 
much  more  convenient  to  work  in  Cartesian  coordinates.  The  expressions 
are  simplified  considerably  by  using  complex  fields  and  complex  variables. 

By  transformation  of  Eqs.  (2.10)  and  (2.11),  we  find  for  the  displace- 
ments. 

ux  - 1uy  = aDjZZ**1  + 02ZX  - (3  - 4v)  D1ZX  (D.l) 

The  stresses  are  of  the  form 

°yy  + ioyx  = '2uX  [Di  ^ + (5i  + D2)zX_1 

(D.2) 

+(A-1)D1  2zx"2]  = *(z,  2) 

oxx  = 2mA  Real  [-2D1  2X_1  + D2  ZX_1  + (A-  1)D]  ZZX'2]  (D.3) 

where  Z = x + iy,  and  the  superior  bar  indicates  the  complex  conjugate,  here 


D1  » C1  + 1C2 
^2  * S * ’^4 


(D.5) 


where  the  C's  are  the  unknown  real  constants  of  Integration,  given  In 
Eqs.  (2.11). 

Referring  to  Eq.  (D.2),  we  see  that  the  condition  of  no  traction  on 
a crack  (x,  yQ  = 0),  x < 0 is  *( x + 1yQ  x - iyg)  = 0,  x < 0. 

In  the  derivation  of  the  perturbed  Williams'  solutions  for  the  wide 
crack,  extensive  use  will  be  made  of  the  derivatives 


(1)"[n  (x  - t)]  {(-l)n  D1Zx'n_1 

3yn  t=0  1 

+ (D,  + D2  - nD1)Zx‘n‘1  + 0] ( x-n-1 )ZZx"n’2} 


(D.6) 


for  the  expansion  of  the  <|>  function  in  a Taylor  series. 

In  addition  to  the  ordinary  Williams'  solutions,  Williams'  logarithmic 
solutions  are  also  used  as  perturbation  functions.  The  n-th  logarithmic 
solution  can  be  written  of  the  form: 


u(n)  _ iu(n)  = nD^2ZX_1  Ln11"1  (Z)  + [xD^Z^1  + D2ZX]  Lnn(Z) 


(D.7) 


- (3  - 4v)D1ZX  Lnn(Z) 


*(n)  - c<y">  ♦ io^  = -2nX  {D1  [n  Lnn-1(Z)  + xLnn(Z)]lX’1 


+ (D]  + 02)  [n  Lnn_1(Z)  + xLnn(Z)]ZX_1 


(D.8) 


+ D1  [n(n-l)  Lnn"2(Z)  + n(2x  - 1)  Lnn_1(Z)  + X2  Lnn(Z)]ZZX-2} 
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ojj*  « 2y  Real  {-201  [n  Lnn_1(2)  + Lnn(2)]ZA_1 

+ 02  [n  Lnn_1(Z)  + XLnn(Z)]ZX_1  (D.9) 

+ [n(n-l)  Lnn'2(Z)  + n(2x  - 1)  Lnn-1(Z)  + X2  Lnn(Z)]Zzx‘2} 


! 
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APPENDIX  E 

GENERAL  FORMS  OF  PERTURBATION  FUNCTIONS  AND 
THEIR  DETERMINATION  FROM  BOUNDARY  CONDITIONS 

E . 1 General  Forms  of  the  Perturbation  Functions 

As  mentioned  in  Section  4.2,  for  the  determination  of  the  perturbed 
solutions  in  the  wide  crack  problem  we  start  from  the  eigenfunctions 
corresponding  to  the  infinitesimal  crack  (zeroth  perturbation),  with 
the  eigenvalues 

x = + ? m = 0,  1,  2,  3 . . . (E.l) 

m - c. 

which  occur  as  double  roots.  The  general  form  of  these  solutions  is 
given  in  polar  coordinates  by  Eqs.  (2.10)  through  Eqs.  (2.13)  and  in 
complex  form  by  Eqs.  (D.l)  through  Eqs.  (D.3).  These  solutions  are 
divided  into  symmetric  and  anti-synmetric  forms.  For  example,  for 
the  half-integer  eigenvalues 

X = + 0.50,  + 1.50,  + 2.50,  + . . . . (E.2) 

m 

The  relations 
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must  hold  in  order  that  the  solutions  be  symmetric  about  the  x-axis.  A 
relation  of  the  form 

S,  • <E-4> 

between  the  other  two  components  of  the  vector  guarantees  homogeneous 
traction  boundary  conditions  along  the  half  line  x < 0,  y a 0,  or  guarantees 
that  we  have  an  eigensolution. 

Similarly  the  eigenvectors  for  the  anti -symmetric  eigenfunctions 
corresponding  to  eigenvalues  given  by  Eq.  (E.2),  are  of  the  form 

Clm  ' C3„  ’ 0 (E'5) 

The  half  line  x < 0,  y = 0 is  freed  of  traction  if 


+ 1)  C 


2m 


(E.6) 


It  can  be  seen  that  the  half  integers  are  double  eigenvalues  for 
the  slit  plane.  The  integers  also  correspond  to  double  eigenvalues  and 
may  be  split  into  symmetric  and  anti-symmetric  solutions.  These  eigen- 
values are  treated  much  the  same  way  as  the  half  integer  ones. 

For  the  x's  given  by  Eq.  (E.2),  one  other  independent  symmetric 
solution  can  be  derived;  this  solution  will,  of  course,  not  be  an  eigen- 
function. A simple  choice  of  the  C's  for  this  solution  is  given  by 
Eq.  (E. 3)  and 


I U. 


1 


C3m  ' 0 


Clm  " 0 


(E.  7) 


which  gives  a (x  < 0,  y = 0)  = 0,  o (x  < 0,  y = 0)  t 0. 

y*  yy 


The  Williams'  logarithmic  solutions  are  other  solution  forms  used  in  the 
perturbation  process.  Here  again,  two  independent  symmetric  solutions  can 
be  derived  by  properly  choosing  the  constants  of  integration.  For  example, 
for  xm  given  by  Eq.  (E.2),  if  the  c's  taken  in  proportion  given  by  Eqs.  (E.3) 
and  (E.4)  the  solution  will  be  symmetric  and  the  tractions  of  the  logarithmic 
term  will  satisfy  the  homogeneous  boundary  condition  along  the  half  line 
x < 0,  y = 0.  If  the  C's  are  chosen  of  the  form  given  by  Eq.  (E.7),  the 
solution  will  give  non-zero  tractions  varying  logarithmically  along  the 
boundary  edges  of  the  corresponding  infinitesimal  crack. 

E . 2 Determination  of  the  Perturbation  Functi ons  from  the  Boundary  Conditions 

The  functions  (Eq.  4.1)  have  to  be  determined  subject  to  the  boundary 
conditions  given  by  Eqs.  (4.4).  As  an  example  let  us  take  a synmetric 
case  and  find  J^for  the  first  perturbation  (k  = 1).  We  have 


* 3 3]*!  + 82*2 


(E.8) 


Here  the  indices  indicating  the  perturbation  number  are  dropped  for  brevity. 
The  p's  are  unknown  coefficients  yet  to  be  determined,  and  the  ^'s  are  the 
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solution  functions  to  be  found  to  satisfy  the  boundary  conditions  given  by 


* 


x < 0 
y = 0 


a*, 


_o 

ay 


x < 0 
y = 0 


(E.9) 


In  Eq.  (E.9)  the  (zeroth  perturbation  function)  ir  known.  It  is 
the  eigenfunction  corresponding  to  the  infinitesimal  crack  with  the 
eigenvalue  given  by  Eq.  (E.2)  and  the  eigenvector  given  by  Eqs.  (E.3) 
and  (E.4).  For  a particular  eigenvalue  say  Xg,  the  right  hand  side  of 
Eq.  (E.9)  can  be  evaluated  using  Eq.  (D.6).  We  obtain 

♦ I » (f0  + ig0)px0'2  (E.  10) 

x < 0 
y = 0 

In  Eq.  (E.10)  one  or  both  of  the  values  fQ  and  g0  may  be  zero 
depending  upon  the  eigenvalue  and  the  symmetric  or  anti -symmetric  nature 
of  the  eigenfunction. 

As  discussed  in  Section  4.2,  one  of  the  functions  in  the  Eq.  (E.9) 
will  be  of  non-eigenfunction  type  with  coefficients  given  by  Eq.  (E.4) 
and  the  other  of  the  form  of  the  Williams'  logarthmic  type  with  the  C’s 
proportional  to  the  eigenvector  of  the  zeroth  eigenfunction.  Evaluating 
these  functions  using  Eq.  (D.4)  and  (D.7)  and  substituting  in  Eq.  (D.8), 
we  obtain 

8^ (fj  + ig-jJp*1"1  + 82^2  + 193)  + (^3  + ■* 93)  Ln(p)]p'2  1 

( E • 1 1 ) 

= (f0  + ig0)pxo‘2 
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The  C's  of  the i^2  function  are  taken  proportional  to  those  of  4>0  which, 
according  to  Section  E.l,  will  give 


f3+ig3  = 0 


(E.12) 


+ ig1)Pxr1  + e2(f2  + 1‘92)pX2"1  = (f0  = i90)pX°"2  (E.13) 

Since  the  solution  function  must  satisfy  the  boundary  conditions  all 
along  the  (x  < 0,  p = |x|)  we  should  choose 


X1 ~ xo'2  » 


X1  = V1 


(E.14) 


X2_  ^ “ xo“2  ’ xo  = xn“ I 


2 A0 


and  finally  Eq.  (E.13)  must  be  true  for  all  values  of  p which  give 


Blfl  + B2f2  = T0 


(E.l 5) 


Blfl  + h*Z  = 9f 


for  the  determination  of  ^ and  ^ 
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APPENDIX  F 

DETERMINATION  OF  DISPLACEMENTS  AND  STRESSES 
FOR  AN  ARBITRARY  POINT  IN  THE  BODY 


F . 1 Formulation  of  the  Problem 

As  mentioned  in  Section  2.1.2,  once  the  unknown  boundary  values  are 
determined,  the  displacements  at  any  point  of  the  boundary  can  be  found 
using  Eq.  (2.5).  The  strain  and  eventually  the  stresses  are  evaluated  by 
differentiation  of  Eq.  (2.4)  or  Eq.  (2.5)  with  respect  to  Xg  and  yg  under 
the  integral  sign.1  We  obtain 


au. . 

[t  (Q)  ( p ,Q ) 


-u.(Q)  -JJ-  (P,Q)]ds 
3x0 

= / Ct.  (Q)  (P.Q) 

Ji  J 3y0 

* 

- ti.(Q)  -U  (P,Q)]ds 


(F.l) 


(F.2) 


The  x0  and  yQ  are  cartesian  coordinate  of  point  P where  the  Kelvin  load 
is  applied.  The  Integration  under  the  Integral  sign  for  the  singular 
integral  is  permissible  [31]. 
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All  the  terms  in  Eqs.  (F.l)  and  (F.2)  are  defined  in  Sections  2.2.1 
and  2.2.2.  The  final  expressions  for  the  Kelvin  kernels  involve  the  chord 
length  r connecting  P and  Q and  Its  derivatives  with  respect  to  x.  y,  n,  s 
coordinates  at  point  Q.  We  have 


r2  ■ (*  - x0)2  * (y  - y0)2 


rn  ' "V  * ”ry 
X*  xrt 


r a fr  mr 
s y x 


y-y0 


(F.3) 


where  «.  and  m are  the  direction  cosine  of  the  normal  to  the  boundary  at 
point  Q.  Furthermore 


ar 


3 rK 

Ln(r)  = — - 

3X0  r 

3rn  rv 

= JL  (-5,  r + mr  ) 

3xft  r y * 


3 rv 

Ln(r)  * — *- 

3y0  r 

(F.4) 

3rn  r* 

- — - (-«  r + mr  ) 

3y0  r y 


a 


^ (t  rx  + mry) 


* 


where 


Rx  - 2rxrn  - 1 


Ry  - fyr„  - ■ 


(F.  7) 


F.2.  Evaluation  of  Integrals  in  the  Singular  Case 

For  the  evaluation  of  the  integrals  of  Eqs.  (2.5),  (F.l)  and  (F.2) 
along  the  intervals  in  which  the  kernels  are  unbounded,  the  density  function 
is  represented  by  a second  order  polynomial  given  by  Eq.  (5.1)  and  the 
kernel  is  expressed  in  terms  of  the  curvature  of  the  boundary  and  its 
derivatives  with  respect  to  the  arc  length  in  the  neighborhood  of  the 
singular  points.  The  expression  for  the  curvature  can  be  written  as 


k"  2 k‘"  3 

k = kQ+k's  + ^ys  + 2j  s + 


(F.8) 


2 

where  k * = — - and  k"  = . Using  the  relation  ~ = k,  we  get 

as  3S*  3S 


= [S  k ds  = 
JO 


lx  I 1 l/»  o 

kgS  + -~2  s + -g-  s + . . . 


(F.9) 


where  the  angle  e is  shown  in  Fig.  13.  Furthermore,  the  chord  length  r and 
its  derivative  with  respect  to  x,  y,  n,  s coordinate  (see  Fig.  13)  can  be 
expanded  using  Eqs.  (F.7)  and  F.8).  We  obtain^ 


1 


See  Forbes  [6]  for  more  detail. 


coSa(S)  . , . . s3  ♦ [ £ . i d _ s«  ♦ (F.,0) 


2 2 


24  8 6 


l • ^2  3 .2,  i « 

sine(s)  = ks  + + [kM  - k ] sq  + 


(F.ll) 


= s - 


ids! . yc  s4  F ^ 


k^  . kkl]$5  + 


120  40  30 


(F.12) 


kid  t + [r  . k3]  i!  . kV  s5  t 
26  24  20 


(F.13) 


r2(s) 


Ln(r) 


= (S2  ♦ 62)  - s2  [1  - idli  . kki£  . (Jil  . id  kkl,  s4] 

12  12  360  45  40  , 


= Ln | s | - {— - + — s3  + . . . . ) 
24  24 


(F.14) 
(F. 15) 


£/%k  o ITI/nIc  ^nk 

-«u  * s[— 2-]  * s2  [-2 2-] 

u 2 8 6 


_ £nk  nukk'  2nk" 

+ s3[-2—  + -9 y_] 


l 

12 

24 

2k- 

m„k4  + 

m0kk" 

,t 5od 

45 

48.8 

48 

72 

(F.16) 


m_  k 0 juk^  mnk 1 
«„  + S[-  -S-]  ♦ s2[ — 2 0_] 

u 2 2 6 


+ sJ[ 

♦ S4[ 


v3 

*0kk" 

V] 

24 

48 

12 

m0k2k ' 

t k4 

iqk 

£0kk" 

. V 

45 

48x8 

48 

72 

(F.17) 


r\>|*r 


( F . 18) 

(F.  19) 

(F.20) 

(F.21) 

(F.22) 


(F.23) 
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. 1 

In  the  singular  Interval,  the  displacements  and  tractions  of  the  Kelvin  ' 

forces  can  be  represented  as  a truncated  laurant  series  with  respect  to  arc 
lengths  using  Eqs.  (F.14)  through  (F.23).  The  general  form  of  the  kernels 
can  be  written  of  the  form  (see  Ref.  6). 

a-l  ? 

g(s)  = “ + «o  + a!s  + v + (f.24) 

i 

Similarly  for  evaluation  of  the  integrals  of  Eqs.  (F.l)  and  (F.2) 
in  the  interval  where  the  kernels  are  unbounded,  the  relations  given  by 
Eqs.  (F.5)  and  ( F . 6)  are  also  expanded  using  Eqs.  (F.14)  through  (F.23), 
the  expansions  are  lengthy  only  few  will  be  given  here 
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] 

I 


au*2  c,  l-t  * 2m  2) 


3X, 


lQV-.  • t.n0  / mnk  _ „ 

1 + C[  "TT  + m0k(2t02  " m02)] 


* *ok2  "’ok‘  , lo3k2  mo3k'  , 2Vo2k*  5mo2lok2 

H 1 ->  ♦ “Z  + — 5 5 + = = — 


~TF~  “6 


t„kk'  m0k"  t03 *kk'  m 3k3  m03k"  m t V 

* C[  * — -j-  * ^ ♦ _a6_ 


1*0  *okk'  *0*0  l*3  t 2 

— nr~  • — j — Is 


it  2 

3t,,(Q)  Aln(B  + mn  ) .r  .t/  u . u \ . _ t./  u u 
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The  other  relations  given  by  Eqs.  (F.5)  and  (F.6)  can  be  evaluated 
in  like  manner.  The  kernels  for  the  strain  equations  can  now  be  written 


in  the  form 


« o a_i  2 
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(F. 29) 


Using  Eq.  (5.6)  the  density  function  f(s)  and  the  kernel,  Eq.  (F.24)  or 
(F.27)  in  Eq.  (5.1)  we  can  easily  evaluate  the  integrals  in  the  singular 
Intervals  for  determination  of  displacements  or  strains. 


PART  1 - Government 


Administrative  and  Liaison  Activities 

Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 
Attn:  Code  474  (2) 

Code  471 
Code  200 

Director 

Office  of  Naval  Research 
Branch  Office 
666  Summer  Street 
Boston,  MA  02210 

Director 

Office  of  Naval  Research 
Branch  Office 
536  South  Clark  Street 
Chicago,  IL  60605 

Director 

Office  of  Naval  Research 
New  York  Area  Office 
715  Broadway  - 5th  Floor 
New  York,  NY  10003 

Director 

Office  of  Naval  Research 
Branch  Office 
1030  East  Green  Street 
Pasadena,  CA  91106 

Naval  Research  Laboratory  (6) 

Code  2627 

Washington,  DC  20375 

Defense  Documentation  Center  (12) 
Cameron  Station 
Alexandria,  VA  22314 


Undersea  Explosion  Research  Division 
Naval  Ship  Research  and  Development 
Center 

Norfolk  Naval  Shipyard 
Portsmouth,  VA  23709 

Attn:  Dr.  E.  Palmer,  Code  177 


Navy  (Con't.) 

Naval  Research  Laboratory 
Washington,  DC  20375 
Attn : Code  8400 
8410 
8430 
8440 
6300 
6390 
6380 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Annapolis,  MD  21402 
Attn : Code  2740 
28 
281 

U.S.  Naval  Weapons  Center 
China  Lake,  CA  93555 
Attn : Code  4062 
4520 

Commanding  Officer 

U.S.  Naval  Civil  Engineering  Laboratory 
Code  L31 

Port  Hueneme,  CA  93041 

Naval  Surface  Weapons  Center 
White  Oak 

Silver  Spring,  MD  20910 
Attn:  Code  R-10 
G-402 
K-82 

Technical  Director 
Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Supervisor  of  Shipbuilding 
U.S.  Navy 

Newport  News,  VA  23607 

U.S.  Navy  Underwater  Sound 
Reference  Division 
Naval  Research  Laboratory 
P.0.  Box  8337 
Orlando,  FL  32806 


(Con' t. ) 


Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  DC  20350 
Attn:  Code  OP-098 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
Washington,  DC  20376 
Attn:  NSP-200 


Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20361 

Attn:  Code  5302  (Aerospace  and  Structures) 
604  (Technical  Library) 

320B  (Structures) 


Commanding  Officer  and  Director 
David  W.  Taylor  Naval  Ship 

Research  and  Development  Center 


Bethesda, 

Attn: 


MD 
Code 


Naval  Air  Development  Center 
Director,  Aerospace  Mechanics 
Warminster,  PA  18974 

U.S.  Naval  Academy 
Engineering  Department 
Annapolis,  MD  21402 


20084 

042 

17 

172 

173 

174 
1800 
1844 
1102.1 

1900 

1901 
1945 
1960 
1962 


Naval  Underwater  Systems  Center 
Newport,  RI  02840 
Attn:  Dr.  R.  Tralnor 


Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22448 
Attn:  Code  G04 
G20 


Naval  Facilities  Engineering  Command 
200  Stovall  Street 
Alexandria,  VA  22332 

Attn:  Code  03  (Research  and  Development)  Technical  Director 


04B 
045 

14114  (Technical  Library) 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 

Attn:  Code  03  (Research  and  Technology) 
037  (Ship  Silencing  Division] 
035  (Mechanics  and  Materials) 

Naval  Ship  Engineering  Center 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  Code  6105G 
6114 
6120D 
6128 
6129 


Mare  Island  Naval  Shipyard 
Vallejo,  CA  94592 

U.S.  Naval  Postgraduate  School 

Library 

Code  0384 

Monterey,  CA  93940 

Webb  Institute  of  Naval  Architecture 
Attn:  Librarian 
Crescent  Beach  Road,  Glen  Cove 
Long  Island,  NY  11542 


Air  Force 


( 


1 

1 

L 

L 

t 


Army 

Commanding  Officer  (2) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709 
Attn:  Mr.  J.  J.  Murray, 
CRO-AA- IP 

Watervllet  Arsenal 
MAGGS  Research  Center 
Watervllet.  NY  12189 

Attn:  Director  of  Research 

U.S.  Army  Materials  and  Mechanics 
Research  Center 
Watertown,  MA  02172 

Attn:  Dr.  R.  Shea,  DRXMR-T 

U.S.  Army  Missile  Research  and 
Development  Center 
Redstone  Scientific  Information 
Center 

Chief,  Document  Section 
Redstone  Arsenal,  AL  35809 

Army  Research  and  Development 
Center 

Fort  Belvolr,  VA  22060 


Commander  WADD 

Wrlght-Patterson  Air  Force  Base 
Dayton,  OH  45433 
Attn:  Code  WWRMDD 

AFFDL  (FDDS) 

Structures  Division 
AFLC  (MCEEA) 

Chief  Applied  Mechanics  Group 
U.S.  Air  Force  Institute  of  Technology 
Wrlght-Patterson  Air  Force  Base 
Dayton,  OH  45433 

Chief,  Civil  Engineering  Branch 
WLRC,  Research  Division 
Air  Force  Weapons  Laboratory 
Klrtland  Air  Force  Base 
Albuquerque,  NM  87117 

Air  Force  Office  of  Scientific  Research 
Bolling  Air  Force  Base 
Washington,  DC  20332 

Attn:  Mechanics  Division 

Department  of  the  Air  Force 
Air  University  Library 
Maxwell  Air  Force  Base 
Montgomery,  AL  36112 


NASA 


National  Aeronautics  and  Space  Administration 

Structures  Research  Division 

Langley  Research  Center 

Langley  Station 

Hampton,  VA  23365 

National  Aeronautics  and  Space  Administration 
Associate  Administrator  for  Advanced 
Research  and  Technology 
Washington,  DC  20546 

Scientific  and  Technical  Information  Facility 
NASA  Representative  (S-AK/DL) 

P.O.  Box  5700 
Bethesda,  MD  20014 


i 


! 

Jt-rj 


Other  Government  Activities 


Commandant 

Chief,  Testing  and  Development  Division 
U.S.  Coast  Guard 
1300  E Street,  NW 
Washington,  DC  20226 

Technical  Director 
Marine  Corps  Development 
and  Education  Conmand 
Quantlco,  VA  22134 

Director  Defense  Research 
and  Engineering 
Technical  Library 
Room  3C128 
The  Pentagon 
Washington,  DC  20301 

Director 

National  Bureau  of  Standards 
Washington,  DC  20034 

Attn:  Mr.  B.  L.  Wilson,  EM  219 

Dr.  M.  Gaus 

National  Science  Foundation 
Environmental  Research  Division 
Washington,  DC  20550 

Library  of  Congress 

Science  and  Technology  Division 

Washington,  DC  20540 

Director 

Defense  Nuclear  Agency 
Washington,  DC  20305 
Attn:  SPSS 

Mr.  Jerome  Persh 
Staff  Specialist  for  Materials 
and  Structures 
OUSDR&E,  The  Pentagon 
Room  3D1089 
Washington,  DC  20301 

Chief,  Airframe  and  Equipment  Branch 
FS-120 

Office  of  Flight  Standards 
Federal  Aviation  Agency 
Washington,  DC  20553 


National  Academy  of  Sciences 
National  Research  Council 
Ship  Hull  Research  Committee 
2101  Constitution  Avenue 
Washington,  DC  20418 
Attn:  Mr.  A.  R.  Lytle 

National  Science  Foundation 
Engineering  Mechanics  Section 
Division  of  Engineering 
Washington,  DC  20550 

Plcatlnny  Arsenal 

Plastics  Technical  Evaluation  Center 
Attn:  Technical  Information  Section 
Dover,  NJ  07801 

Maritime  Administration 
Office  of  Maritime  Technology 
14th  and  Constitution  Ave.,  NW 
Washington,  DC  20230 

Maritime  Administration 
Office  of  Ship  Construction 
14th  and  Constitution  Ave.,  NW 
Washington,  DC  20230 


J 


PART  2 - Contractors  and  Other  Technical  Collaborators 


Universities 


Dr.  J.  Tinsley  Oden 
University  of  Texas  at  Austin 
345  Engineering  Science  Building 
Austin,  TX  78712 

Professor  Julius  Mlklowltz 
California  Institute  of  Technology 
Division  of  Engineering 
and  Applied  Sciences 
Pasadena,  CA  91109 

Dr.  Harold  Llebowltz,  Dean 
School  of  Engineering  and 
Applied  Science 
George  Washington  University 

Professor  Eli  Sternberg 
California  Institute  of  Technology 
Division  of  Engineering  and 
Applied  Sciences 
Pasadena,  CA  91109 

Professor  Paul  M.  Naghdl 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  CA  94720 


Professor  J.  Kempner 
Polytechnic  Institute  of  New  York 
Department  of  Aerospace  Engineering  and 
Applied  Mechanics 
333  Jay  Street 
Brooklyn,  NY  11201 

Professor  J.  Klosner 
Polytechnic  Institute  of  New  York 
Department  of  Aerospace  Engineering  and 
Applied  Mechanics 
333  Jay  Street 
Brooklyn,  NY  11201 

Professor  R.  A.  Schapery 
Texas  A&M  University 
Department  of  Civil  Engineering 
College  Station,  TX  77843 

Professor  Walter  D.  Pllkey 
University  of  Virginia 
Research  Laboratories  for  the 
Engineering  Sciences 
School  of  Engineering  and 
Applied  Sciences 
Charlottesville,  VA  22901 


Professor  A.  J.  Durelll 
Oakland  University 
School  of  Engineering 
Rochester,  MI  48063 

Professor  F.  L.  DiMagglo 
Columbia  University 
Department  of  Civil  Engineering 
New  York,  NY  10027 

Professor  Norman  Jones 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge,  MA  02139 

Professor  E.  J.  Skudrzyk 
Pennsylvania  State  University 
Applied  Research  Laboratory 
Deoartment  of  Physics 
State  College,  PA  16801 


Professor  K.  D.  Willmert 
Clarkson  College  of  Technology 
Department  of  Mechanical  Engineering 
Potsdam,  NY  13676 

Dr.  Walter  E.  Haisler 
Texas  A&M  University 
Aerospace  Engineering  Department 
College  Station,  TX  77843 

Dr.  Hussein  A.  Kamel 
University  of  Arizona 
Department  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  A Z 85721 

Dr.  S.  J.  Fenves 
Carnegle-Mellon  University 
Department  of  Civil  Engineering 
Schenley  Park 
Pittsburgh,  PA  15213 


Universities  (Con't.) 


Dr.  Ronald  L.  Huston 
Department  of  Engineering  Analysis 
University  of  Cincinnati 
Cincinnati,  OH  45221 

Professor  6.  C.  M.  SI h 
Lehigh  University 
Institute  of  Fracture  and 
Solid  Mechanics 
Bethlehem,  PA  18015 

Professor  Albert  S.  Kobayashl 
University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  WA  98105 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
State  University 

Department  of  Engineering  Mechanics 
Blacksburg,  VA  24061 

Professor  A.  C.  Erlngen 
Princeton  University 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  NJ  08540 

Professor  E.  H.  Lee 

Stanford  University 

Division  of  Engineering  Mechanics 

Stanford,  CA  94305 

Professor  Albert  I.  King 
Wayne  State  University 
Biomechanics  Research  Center 
Detroit,  MI  48202 

Dr.  V.  R.  Hodgson 
Wayne  State  University 
School  of  Medicine 
Detroit,  MI  48202 

Dean  B.  A.  Boley 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  IL  60201 


Professor  P.  6.  Hodge,  Jr. 
University  of  Minnesota 
Department  of  Aerospace  Engineering 
and  Mechanics 
Minneapolis,  MN  55455 

Dr.  D.  C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbana,  IL  61801 

Professor  N.  M.  Newmark 
University  of  Illinois 
Department  of  Civil  Engineering 
Urbana,  IL  61803 

Professor  E.  Relssner 
University  of  California,  San  Diego 
Department  of  Applied  Mechanics 
La  Jolla,  CA  92037 

Professor  William  A.  Nash 
University  of  Massachusetts 
Department  of  Mechanics  and 
Aerospace  Engineering 
Amherst,  MA  01002 

Professor  G.  Herrmann 
Stanford  University 
Department  of  Applied  Mechanics 
Stanford,  CA  94305 

Professor  J.  D.  Achenbach 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  IL  60201 

Professor  S.  B.  Dong 
University  of  California 
Department  of  Mechanics 
Los  Angeles,  CA  90024 

Professor  Burt  Paul 
University  of  Pennsylvania 
Towne  School  of  Civil  and 
Mechanical  Engineering 
Philadelphia.  PA  19104 


P 


Uni versl ties  (Con't.) 


Professor  H.  W.  Liu  Dr.  Samuel  B.  Batdorf 

Syracuse  University  University  of  California 

Department  of  Chemical  Engineering  School  of  Engineering 
and  Metallurgy  and  Applied  Science 

Syracuse,  NY  13210  Los  Angeles,  CA  90024 


Professor  S.  Bodner 
Technlon  R&D  Foundation 
Haifa,  Israel 

Professor  Werner  Goldsmith 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  CA  94720 

Professor  R.  S.  Rlvlin 
Lehigh  University 
Center  for  the  AodII cation 
of  Mathematics 
Bethlehem,  PA  18015 

Professor  F.  A.  Cozzarelll 
State  University  of  New  York  at  Buffalo 
Division  of  Interdlsclpl Inary  Studies 
Karr  Parker  Engineering  Building 
Chemistry  Road 
Buffalo,  NY  14214 

Professor  Joseph  L.  Rose 
Drexel  University 

Department  of  Mechanical  Engineering 
and  Mechanics 
Philadelphia,  PA  19104 

Professor  B.  K.  Donaldson 
University  of  Maryland 
Aerospace  Engineering  Department 
College  Park,  MD  20742 

Professor  Joseph  A.  Clark 
Catholic  University  of  America 
Department  of  Mechanical  Engineering 
Washington,  DC  20064 

Professor  T.  C.  Huang 
University  of  Wlsconsln-Madlson 
Department  of  Engineering  Mechanics 
Madison,  WI  53706 


Professor  Isaac  Fried 
Boston  University 
Department  of  Mathematics 
Boston,  MA  02215 

Professor  Michael  Pappas 
New  Jersey  Institute  of  Technology 
Newark  College  of  Engineering 
323  High  Street 
Newark,  NJ  07102 

Professor  E.  Krempl 
Rensselaer  Polytechnic  Institute 
Division  of  Engineering 
Engineering  Mechanics 
Troy,  NY  12181 

Dr.  Jack  R.  Vinson 
University  of  Delaware 
Department  of  Mechanical  and  Aerospace 
Engineering  and  the  Center  for 
Composite  Materials 
Newark,  DL  19711 

Dr.  Dennis  A.  Nagy 
Princeton  University 

School  of  Engineering  and  Applied  Science 
Department  of  Civil  Engineering 
Princeton,  NJ  08540 

Dr.  J.  Duffy 
Brown  University 
Division  of  Engineering 
Providence,  RI  02912 

Dr.  J.  L.  Swedlow 
Carnegle-Mellon  University 
Department  of  Mechanical  Engineering 
Pittsburgh,  PA  15213 

Dr.  V.  K.  Varadan 

Ohio  State  University  Research  Foundation 
Department  of  Engineering  Mechanics 
Columbus,  OH  43210 


Universities  ( Con ' t . ) 


Industry  and  Research  Institutes  (Con't.) 


Or.  Z.  Hashln 

University  of  Pennsylvania 
Department  of  Metallurgy  and 
Materials  Science 
College  of  Engineering  and 
Applied  Science 
Philadelphia,  PA  19104 

Dr.  Jackson  C.  S.  Yang 
University  of  Maryland 
Department  of  Mechanical  Engineering 
College  Park,  MD  20742 

Professor  T.  Y.  Chang 
University  of  Akron 
Department  of  Civil  Engineering 
Akron,  OH  44325 

Professor  Charles  W.  Bert 
University  of  Oklahoma 
School  of  Aerospace,  Mechanical, 
and  Nuclear  Engineering 
Norman,  OK  73019 

Professor  Satya  N.  Atlurl 
Georgia  Institute  of  Technology 
School  of  Engineering  Science  and 
Mechanics 

Atlanta.  GA  30332 

Professor  Graham  F.  Carey 
University  of  Texas  at  Austin 
Department  of  Aerospace  Engineering 
and  Engineering  Mechanics 
Austin,  TX  78712 

Industry  and  Research  Institutes 

Dr.  Jackson  C.  S.  Yang 
Advanced  Technology  and  Research,  Inc. 
10006  Green  Forest  Drive 
Adel  phi,  MD  20783 

Dr.  Norman  Hobbs 
Kaman  AvIDyne 
Division  of  Kaman 
Sciences  Corp. 

Burlington,  MA  01803 


Argonne  National  Laboratory 
library  Services  Department 
9700  South  Cass  Avenue 
Argonne,  IL  60440 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
1033  Massachusetts  Avenue 
Cambridge,  MA  02138 

Dr.  V.  Godlno 

General  Dynamics  Corporation 
Electric  Boat  Division 
Groton , CT  06340 

Or.  J.  E.  Greenspon 

J.  G.  Engineering  Research  Associates 

3831  Menlo  Drive 

Baltimore,  MD  21215 

Dr.  K.  C.  Park 

Lockheed  Missile  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

Newport  News  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  VA  23607 

Dr.  W.  F.  Bozich 
McDonnell  Douglas  Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  TX  78284 

Dr.  R.  C.  DeHart 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  TX  78284 

Dr.  M.  L.  Baron 
Weidllnger  Associates 
110  East  59th  Street 
New  York,  NY  10022 


Industry  and  Research  Institutes  (Con't.) 


Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

Mr.  William  Caywood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  MD  20810 

Or.  Robert  E.  Nlckell 
Pacifica  Technology 
P.0.  Box  148 
Del  Mar,  CA  92014 


Dr.  M.  F.  Kannlnen 

Battel le  Columbus  Laboratories 

505  King  Avenue 

Columbus,  OH  43201 

Dr.  G.  T.  Hahn 

Battel le  Columbus  Laboratories 
505  King  Avenue 
Columbus,  OH  43201 

Dr.  A.  A.  Hochreln 
Daedalean  Associates,  Inc. 
Sprlnglake  Research  Center 
15110  Frederick  Road 
Woodbine,  MD  21797 

Mr.  Richard  Y.  Dow 
National  Academy  of  Sciences 
2101  Constitution  Avenue 
Washington,  DC  20418 

Mr.  H.  L.  Kington 
A1 research  Manufacturing  Company 
of  Arizona 
P.0.  Box  5217 
111  South  34th  Street 
Phoenix,  AZ  85010 


Dr.  M.  H.  Rice 

Systems,  Science,  and  Software 
P.0.  Box  1620 
La  Jolla,  CA  92037 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RAO*  Of  Knfrwd) 


READ  INSTRUCTIONS 


REPORT  DOCUMENT ATION  PAGE 


PEKOHE  COMPLETING  RORM 
J RECIPIENT'S  CATALOG  NUMIEN 


TYPE  QF  Rl 

^Technica 


Boundary  Integral  Equation  Solution  of 
Plane  Elasticity  Problems  with  High 
Stress  Concentrations  • \ 


CONTRACT  OR  GRANT  NUMBER^*) 


7 AUTHOR^ »)_ 


Hassan/ii 
Arthur* R 


100014- 75-C-0164 


Robinson 


10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 


(FORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  Civil  Engineering 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  Illinois  61801 


Project  No.  NR  064-183 


II  CONTROLLING  OFFICE  NAME  ANO  ADDRESS  /^' 

Office  of  Naval  Research,  CODE  N00014  ( I 

Department  of  the  Navy  V_ 

Arlington,  Virginia  22217 

n MONITORING  AGENCY  NAME  A ADDRESS^**  dlllaront  feom  Controlling  Oil lea) 

Material  Science  Division 

Structural  Mechanics  Program  (CODE  474) 

Office  of  Naval  Research  (800  Quincy  Street) 

22217 


IS  SECURITY  CLASS  (of  thla  report) 

UNCLASSIFIED 


1S«  DC  CL  ASSIFIC  ATION  DOWNGRADING 
SCHEDULE 


Arlington,  Virginia 


IE  DISTRIBUTION  STATEMENT  (ol  thl,  N»porl) 


Approved  for  public  release:  Distribution  Unlimited 


17  DISTRIBUTION  STATEMENT  (ol  Iho  abstract  antarad  In  Mock  20,  It  dltlacant  from  Rapor(\ 


ft  KEY  WOROS  ( Continue  on  rararaa  ml  da  If  nacaaamry  and  Identify  by  block  number) 

Plane  Elasticity 

Boundary  Integral  Equations 

Stress  Concentrations 

Notches 

Fillets 

A numerical  method  for  determination  of  stresses  in  two-dimensional 
elastic  bodies  with  high  stress  concentrations  is  presented.  Major 
emphasis  Is  placed  on  bodies  with  a notch  having  a fillet  of  small  radius 
and  bodies  with  a crack  of  small  width.  These  static  boundary  value 
problems  are  formulated  in  terms  of  boundary  integral  equations  of  a type 
used  previously  by  Barone  and  Robinson  for  sharp  notches  and  cracks.  . 


security  classification  of  this  page  cBB«n  d»i« 


f cuwty  clarification  or  thu  p»o«(wm»  d««  *»*»q 


20^. (Continued) 

For  the  fillet  problem  a small  Inner  region  containing  the  fillet  and 
bounded  by  a circle  Is  separated  out  and  analyzed  numerically  under  the 
loading  system  of  each  of  the  Williams'  solution  of  the  corresponding 
sharp  notch.  The  results  are  then  used  to  develop  analytical  solutions  for 
the  Intermediate  region  adjacent  to  the  fillet  region.  In  this  way  the 
details  of  the  boundary  configuration  of  the  fillet  Is  reflected  in  a set 
of  generalized  displacements  which  characterize  the  Intermediate  field. 


For  the  solution  of  the  whole  body  Including  the  loaded  boundary, 
kernels  developed  by  Barone  and  Robinson  are  used  to  pick  out  the 
generalized  displacements  of  the  Intermediate  region. 


The  solution  of  the  wide  crack  problem  is  similar  except  that  a 
perturbation  scheme  Is  required  to  ensure  homogeneous  boundary  conditions 
at  the  actual  edges  of  the  crack. 
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